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Abstract 

Percolation is the simplest fundamental model in statistical mechanics that exhibits phase 
transitions signaled by the emergence of a giant connected component. Despite its very 
simple rules, percolation theory has successfully been applied to describe a large variety 
of natural, technological and social systems. Percolation models serve as important 
universality classes in critical phenomena characterized by a set of critical exponents 
which correspond to a rich fractal and scaling structure of their geometric features. We 
will hrst outline the basic features of the ordinary model. 

Over the years a variety of percolation models has been introduced some of which 
with completely different scaling and universal properties from the original model with 
either continuous or discontinuous transitions depending on the control parameter, di¬ 
mensionality and the type of the underlying rules and networks. We will try to take a 
glimpse at a number of selective variations including Achlioptas process, half-restricted 
process and spanning cluster-avoiding process as examples of the so-called explosive per¬ 
colation. We will also introduce non-self-averaging percolation and discuss correlated 
percolation and bootstrap percolation with special emphasis on their recent progress. 
Directed percolation process will be also discussed as a prototype of systems displaying 
a nonequilibrium phase transition into an absorbing state. 

In the past decade, after the invention of stochastic Lowner evolution (SLE) by Oded 
Schramm, two-dimensional (2D) percolation has become a central problem in probability 
theory leading to the two recent Fields medals. After a short review on SLE, we will 
provide an overview on existence of the scaling limit and conformal invariance of the 
critical percolation. We will also establish a connection with the magnetic models based 
on the percolation properties of the Fortuin-Kasteleyn and geometric spin clusters. As 
an application we will discuss how percolation theory leads to the reduction of the 3D 
criticality in a 3D Ising model to a 2D critical behavior. 

Another recent application is to apply percolation theory to study the properties 
of natural and artificial landscapes. We will review the statistical properties of the 
coastlines and watersheds and their relations with percolation. Their fractal structure 
and compatibility with the theory of SLE will also be discussed. The present mean 
sea level on Earth will be shown to coincide with the critical threshold in a percolation 
description of the global topography. 
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1. Introduction 

Consider a simple electric circuit consisting of a voltage source and a bulb connected to 
an insulating hexagonal honeycomb lattice of size N. Now imagine we start occupying the 
lattice in a random way by n number of metallic hexagonal plaquettes. For small values 
of fraction of occupation p = n/N, lower than a critical threshold Pc, there appear some 
small metallic clusters but the two opposite sides of the lattice still remain disconnected 
leading to the lack of charge flow in the circuit and thus the bulb remains off—see 
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Fig. [3 A metallic cluster is defined as a set of occupied sites that can be traversed 
by jumping from neighbor to occupied neighbor. In Fig. [3 the different colors simply 
label the different isolated clusters and have no other significance. By increasing the 
number of the randomly distributed metallic plaquettes (and increasing the occupation 
probability p indeed), the average size of the metallic clusters increases. Once it exceeds a 
certain threshold value Pc, a spanning cluster emerges that closes the circuit and the bulb 
suddenly lights up. This transition from insulator to the metallic phase in two dimensions 
(2D) exemplifies one of the simplest and fundamental classes of phase transitions in 
statistical physics called ’’percolation,” which can generally be defined in any dimension 
d. 

The idea of percolation model was first effectively considered by Flory in 1940s mw- 
However, the study of the model as a mathematical theory, dates back to 1954 [3], when 
engineer Simon Broadbent and mathematician John Hammersley, one concerned with the 
design of carbon filters for gas masks, put their heads together to deal with “The stepwise 
spreading of a fluid or individual particles through a medium following a random path 
in which each link is either open or shut, according to a specihed statistical proportion. 
The spreading process is therefore arrested at many sites...” [5]. A long path along which 
spreading is not interrupted constitutes an infinite (or spanning) cluster. Broadbent and 
Hammersley proposed the concept of a percolation threshold above which the links form 
an infinite cluster with high probability. 

Percolation theory was then popularized in the physics community and intensively 
studied by physicists [giTiiiiigiiioiiiiiiiaiiiiiniiiiiiiiiiiii. It has been found to have 
a broad applications to diverse problems as understanding conducting materials [HIS], 
the fractality of coastlines HOllIIl, networks IS1I1SII21I, turbulence [H US] , magnetic 
models [571 HSl IS] 1301131] , colloids [H |33] , growth models [33] , retention capacity and 
watersheds of landscapes [33 1301133, the spin quantum Hall transition [J8j and SU(3) 
lattice gauge theory |89] . From a mathematical point of view percolation is also attractive 
because it exhibits relations between probabilistic and algebraic/topological properties 
of graphs. Although a significant amount of research has been done in the field there 
still exist many unsolved problems m- 

One of the active areas in percolation research which has still open problems is to 
find the percolation thresholds Pc, as a fundamental characteristic of percolation theory, 
both exactly and by simulation. The values of percolation thresholds are not universal 
and generally depend on the structure of the lattice and dimensionality, and are believed 
to achieve their mean-field values only in the limit of infinite dimension [3- Finding 
mathematical rigorous proofs for thresholds and bounds has also been of research interest 
in [4n l42l l43] . Exact thresholds in 2D for the square, triangular, honeycomb and related 
lattices were found using the star-triangle transformation |44j . It has been shown in 
[35] that thresholds can be found for any lattice that can be represented as a self-dual 
3-hypergraph (that is, decomposed into triangles that form a self-dual arrangement). 
It is also shown in [35] that for any lattice which can geometrically be represented as 
an isoradial graph, it is possible to find the thresholds. This leads to a wide class of 
exact thresholds and provides a proof [33 of the conjecture [48] for the threshold of the 
checkerboard lattice. However, the exact value of thresholds for many systems of long 
interest are still missing [30] . 

For our mentioned example at the beginning of the section, it is equivalent to a 2D 
site percolation problem on a triangular lattice whose sites are placed at the center of the 

2 




Figure 1: Schematic illustration of a simple circuit including a bulb and a voltage source connected to 
an insulator honeycomb lattice which can be furnished by a fraction p of metallic plaquettes. Once p 
reaches a critical threshold pc, a spanning metallic cluster emerges which connects the two sides of the 
lattice and the bulb suddenly lights up. 


plaquettes with six number of nearest neighbors. The value of the threshold for this case 
is exactly known to be Pc = 1/2 in the infinite system size limit. In fact, if our system 
size was infinitely large or equivalently, we could have infinitely small plaquettes, then 
we would need N —> oo number of plaquettes to totally furnish the lattice. In that case, 
the percolation probability defined as the probability to have an infinite cluster Q would 
be a step function around Pc which takes the value 1 for p > p^ and zero if p < Pc—see 
Fig-i Pa p = Pc, it is believed that with probability one there is no infinite cluster, but 
there are typically some very large clusters nearby [15] . 

One may notice that in our example above, we could instead have occupied the edges of 
the sites (or the bonds) with some metallic rods rather than occupying the plaquettes 
themselves. This would then change the problem to one that is known as bond percolation 
on a 2D honeycomb lattice with three number of nearest neighbors. This approach 


^However, one should be careful with the definition of ’’infinite”. 
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changes the percolation threshold to 1 — 2sin(7r/18) « 0.65271 [33], but does not affect 
the other fundamental properties, we shall return to this point later. 

The other point is that in reality the insulator lattice in our experiment is not infinite of 
course, and it has a finite size. The problem of how to deal with finite size lattices is known 
as finite size scaling. It is a useful introduction to the style of theoretical argument that 
is often used in percolation theory [16j . In such situations if we repeat our experiment 
many times for a given occupancy p with different realizations of randomness to estimate 
the percolation probability, we would find a smooth function as represented in Fig. 
rather than a step function around Pc- This means that in reality we can get connectivity 
even at very much less than the percolation threshold or not get it even at a much higher 
occupancy. As the size of the system gets larger the scatter around the sharpness would 
reduce until we return to the plot for the infinite system size. 



Figure 2: Schematic plot of the percolation probability as a function of the occupancy. As the system 
size N goes to infinity, Ps tends to the step function around the critical threshold pc. For very special 
systems all graphs for finite N cross at a single point, but in general, because of finite-size corrections, 
the crossing does not occur at a single point for finite systems |5Q1151| . 


The sudden onset of a spanning cluster at a particular value of the occupation prob¬ 
ability along with a number of characteristic features make the percolation transition 
a nontrivial critical behavior. This criticality belongs to a large family of critical phe¬ 
nomena with common remarkable features and forms an important universality class 
characterized by a number of scaling laws and critical exponents. The notion of uni¬ 
versality means that the large scale behavior of critical systems can be described by 
relatively simple mathematical relationships which are entirely independent of the small 
scale construction. This property enables one to study and understand the behavior 
of a very wide range of systems without needing to know much about the details. For 
example, the universal properties of the percolation model are entirely independent of 
the type of lattice (e.g., hexagonal, triangular or square, etc.) or whether it is site or 
bond percolation; they only depend on the dimensionality of the system. 

Another example of a prototype model in critical phenomena is the Ising model as a 
mathematical model of a magnet which undergoes a phase transition between a ferro¬ 
magnetic ordered phase and a paramagnetic disordered phase at a particular temperature 
Tc, the Curie temperature, in more than one dimension. The net magnetization which is 
the first derivative of the free energy with respect to the applied magnetic field strength, 
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is the order parameter which distinguishes between these two phases i.e., it takes the 
value 0 in the paramagnetic phase and increases continuously from zero as the tempera¬ 
ture is lowered below the Curie point. Such transitions in which the order parameter is 
continuous but second derivative of the free energy (like magnetic susceptibility and heat 
capacity in the Ising model) exhibits a discontinuity are called second-order or continuous 
phase transitions. First-order (discontinuous) transitions involve a discontinuous change 
in their order parameter. Both percolation and Ising models display a continuous phase 
transition and as will be made clear later in Sec. they are closely related to each other 
I13I1H1- Therefore we have to be able to define an order parameter for the percolation 
model as in the Ising model. It seems not to be so difficult: not all occupied sites are in 
the infinite (or spanning) cluster, thus if we look at the probability Poo(p) that an occu¬ 
pied site is in the infinite cluster for a given occupancy p, then this must be zero (since 
there is no spanning cluster) below the percolation threshold and increases continuously 
^ut of course with singular derivatives) as one enters the supercritical (connected) phase 
[j —see Fig. Continuity of the strength probability Poo(p), as the main macroscopic 
observable in percolation, at pc is an open mathematical problem in the general case, 
but it is known to hold rigorously in 2D and d > 19 [52] using lace expansion methods. 
The conjecture that Poo(p = Pc) = 0 for 3 < d < 18 remains however one of the open 
problems in the field [55] . 



Figure 3: Schematic diagram of the percolation strength Poo(p) and the correlation length g(p) as a 
function of the occupancy p for a system of infinite size. 

Percolation is a simple model of robustness and stability exhibiting a robust con¬ 
tinuous transition in all dimensions. It accounts as a fundamental step in dealing with 
more complex models and even dynamical processes occurring on the networks. There 
have been introduced various modifications of the ordinary percolation to model dif¬ 
ferent statistical phenomena or even at the level of mathematical curiosity in building 
up a modified version with a possibly discontinuous phase transition |54l ES] EU l57] . 
Explosive percolation [55], for example, is one of the most challenging cases with a 
seemingly mild modification of standard percolation model which was first, surprisingly, 
claimed to exhibit a discontinuous phase transition, in contrast to the ordinary perco¬ 
lation. It then has been realized that the finite size effects make it difficult to numeri- 


^The system is called to be in the subcritical (resp. supercritical) phase if p < pc (resp. p > pc). 
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cally realize if the explosive percolation in random graphs is continuous or discontinuous 
[551 [501 ion [001 [001 [oll [051 [ok [071 [051 [001 [701 [7T] . Finally, after a numerical observation 
|62| , there has been given a mathematical proof m in favor of the continuous phase tran¬ 
sition. The type of transition for the explosive percolation in Euclidean space joniiTO] 
has not yet been clarified. 

Other variants and modihcations of percolation models will be discussed in Sec. 
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Figure 4: (a) Illustration of the transformed image under the special conformal mapping f{z) = zl{2 — z). 
Left panels in (b) and (c) show the homogeneous spin configurations of a 2D Ising model on a square 
lattice above and at the Curie point, respectively. Under the conformal transformation f{z), the spin 
configuration above Tc is no longer homogeneous while the one exactly at Tc still is and looks statistically 
like the original configuration due to its conformal invariant symmetry. 


Looking back at Figure we see that all graphs of different size cross at a single 
point which is the percolation threshold. This means that at an exactly certain critical 
occupancy the percolation probability becomes scale invariant. Scale invariance means 
that if we zoom in or out and look at the system with an arbitrary window size, then the 
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picture still has the same statistical properties with the same resulting physics. This is 
intimately related to one of the most important methods developed in theoretical physics 
called renormalization group, which studies the behavior of a system under scale trans¬ 
formations m- Renormalization group has been very successful to provide very good 
approximate (or sometimes exact) values of the critical exponents and thresholds for the 
percolation models. 

In two dimensions it turns out that an even stronger symmetry holds at the critical point: 
if we blow up different parts of a figure of percolation clusters by different magnification 
factors (as long as angles are preserved) then statistically the picture once again looks 
the same. This property is called conformal invariance, and with this assumption, the¬ 
oreticians have been able to derive many important properties of critical systems in the 
past few decades [13 Eg. In fact the heuristics behind conformal invariance is a natural 
generalization of rotation and scale invariance. Conformal invariance is much more pow¬ 
erful in two dimensions where it is associated with the theory of analytic functions of a 
complex variable. 

Figure [^visualizes how such transformations work. Consider a special conformal trans¬ 
formation given by f{z) = zl{2 — z) which acts in a way that is shown in Fig. 13 a) on 
the regions of different color in the complex z-plane, the map of each domain of specific 
color in the z-plane at left is shown with the same corresponding color at right. If we 
have an Ising model on a square lattice and color the spin up (down) clusters in black 
(white), then we see that the figures both above and at the Curie point seem to be 
homogeneous—see left panels in Fig. j^b) andj^c). However when we look at the trans¬ 
formed spin configurations under the conformal mapping f{z), we see that for T > the 
figure is no longer homogeneous, while the one at the critical point still is, and looks sta¬ 
tistically the same as the original spin configuration at due to its conformal invariant 
symmetry. 

The emergence of conformal symmetry at the critical point is however more myste¬ 
rious. This seems to be a generic feature of criticality but why this happens is not fully 
understood m- Motivated by numerical experiments [771 EH) in which concluded that 
crossing probabilities should have a universal scaling limit, which is conformally invari¬ 
ant (a conjecture that is attributed to Michael Aizenman [IH]), John Cardy used these 
heuristic ideas in 1992 to give an explicit formula that determines the exact values of the 
crossing probabilities between the opposite sides of a conformal rectangle filled with a 
conformally invariant infinitesimal lattice m- More recently Smirnov rigorously proved 
[50] that the Cardy’s conjecture holds for the continuum limit of site percolation on a 
triangular lattice. But how can one characterize the continuum limit of a lattice model? 

The continuum limit of a lattice model is often difhcult to be captured mathematically. 
In this limit the lattice spacing a is sent to zero, and new sites are constantly added to the 
lattice during this contraction so that the lattice does not eventually vanish but continues 
to fill the original domain that it occupied |81j . The continuum limits of most lattice 
models are believed to converge to the quantum field theories (QFTs). At the critical 
point, since the correlation length ^ is infinite (^ S> a)(Fig. the lattice models must 
be invariant under scale transformations. This property along with the invariance under 
translations and rotations imply (under broad conditions [76| 1 the conformal invariance, 
and indeed, suggest that the continuum limit of critical lattice models should be given by 
conformal field theories (CFTs). Furthermore, only certain CFTs, usually the minimal 
models, have been observed to possess the right structure to describe a critical lattice 
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model in two dimensions. Due to the relatively few number of such theories, models with 
the same macroscopic but different microscopic properties are presumed to have identical 
continuum limits which correspond to the same CFT characterized by the value of the 
central charge c. This is a restatement of the notion of universality discussed earlier. The 
critical percolation and Ising models are famous examples of CFTs with central charge 
c = 0 and c = 1 / 2 , respectively. 

A relatively new method to describe the continuum limit of the critical lattice models 
is called stochastic Lowner evolution (SLE^) invented by Oded Schramm around 2000 
[55] — to review SEE see |53| and |51|. Six years later, in 2006, Wendelin Werner received 
the Fields medal for his contributions to the development of SEE and related subjects. 
Theory of SEE^ is a subject of probability theory that generates planar random curves 
with conformally invariant probability measures in a domain with a boundary. The 
diffusivity 0 < k < 8 is a real parameter that classifies different conformally invariant 
interfaces. For example, the scaling limit of a percolation cluster boundary (hull) is 
proven by Smirnov |5D] to be given by SEEg. It is also shown in [5S] that the interfaces 
in the planar critical Ising model and its random-cluster representation converge to SEE 3 
and SEEig/ 3 , respectively. In 2010, Smirnov was awarded the Eields medal for the proof 
of conformal invariance of percolation and the planar Ising model in statistical physics. 
SEE has soon found many applications and turned out to describe the vorticity lines in 
turbulence |S5| , domain walls of spin glasses [saiHiisi], the nodal lines of random 
wave functions |90| [9T] , the iso-height lines of random grown surfaces [M] ES] [Ml EH EH 
EH, the avalanche lines in sandpile models [98] and the coastlines and watersheds on Earth 
Em [Tun[TU5| . Among which, SEE could provide quite unexpected connections between 
some features of interacting systems and ordinary uncorrelated percolation [HES]. 

In 1969, Fortuin and Kasteleyn (FK) [^ [55] 110311104j found an interesting mapping 
between the g—state Potts model, which includes the Ising model for q = 2, and a corre¬ 
lated bond-percolation model called the random-cluster model. This yielded a geometric 
representation of the partition function for Potts models in terms of the statistics of 
the random clusters. The uncorrelated bond percolation model itself can be recovered 
by taking the limit g —>■ 1 in the FK formalism. This representation has also become 
a key point to derive many analytical results in percolation m- Swendsen and Wang 
una, and then Wolff [106j , have exploited this mapping to devise extraordinarily efficient 
Monte Carlo algorithms based on nonlocal cluster update for Potts models having far 
less critical slowing down than the standard single-spin-update algorithms. It turned out 
that for those values of q for which the model undergoes a continuous phase transition, 
the percolation of FK clusters occurs exactly at the critical temperature and their fractal 
structure encodes the complete critical behavior. It can be shown that there is a one- 
to-one correspondence between different thermodynamic quantities and their geometric 
counterparts based on the statistical and fractal properties of FK clusters. 

A major breakthrough in statistical physics was the exact solution of the Ising model 
in two dimensions [107] . Onsager gave in 1944 a complete solution of the problem in zero 
external magnetic field. But in three dimensions, Istrail has shown |1()8| that essentially 
all versions of the Ising model are computationally intractable across lattices and thus 
the 3D Ising model, in its full, is NP-complete. We will show that is possible to map, 
at least, the criticality of a 3D Ising model onto a 2D cross section of the model. This 
mapping provides a dimensional reduction in the geometrical interpretation of the 3D 
Ising model. Eoosely speaking, for the Ising model on square lattice, there exists an 
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alternative description of the partition function as a sum over all curves surrounding 
the geometric spin clusters (rather than the FK clusters) weighted by their length. The 
continuum limit of the model in this case is shown to be well-defined and is given by the 
theory of free Majorana fermions ung. Thus our attempt to recast the 2D Ising model 
as a theory of immersed curves seems successful. But for an Ising model on a cubic 
lattice, the boundaries between geometric spin clusters are closed surfaces, not curves. 
It is, however, by no means clear how to take the continuum limit of this lattice surface 
theory pg HTMITTOlfTmiTT^ . What we will suggest in Sec. is to replace the lattice 
surface theory in 3D criticality, by a theory of immersed curves on a 2D cross section of 
the original lattice. 

Percolation theory also provides a suitable platform to study the properties of real and 
artificial landscapes. A landscape is a height prohle usually defined on a square lattice 
where each cell’s elevation value at position x represents the average elevation over the 
entire footprint of the cell (site). Now imagine that the water is dripping uniformly over 
the landscape and fills it from the valleys to the mountains, letting the water flow out 
through the open boundaries. During the raining, watershed lines may also form which 
divide the landscape into different drainage basins—see also Fig. Watersheds play a 
fundamental role in geomorphology in e.g., water management [I13j and landslide and 
flood prevention m- For a given landscape represented as a digital elevation map 
(DEM), it is possible to determine the watershed lines based on the iterative application 
of invasion percolation |115j . It is found |I16j that the main watershed line generated 
on random uncorrelated landscapes are self-similar with fractal dimension d/ = 1.2168 ± 
0.0005. It has been also shown that a watershed lines is statistically compatible with the 



Figure 5: Fractal structure of a natural watershed at the north of Tehran, Darabad. The main watershed 
is the (blue) curve connecting the two points marked by ★. 


family of conformally invariant SLE„ curves with k = 1.734±0.005 [93]. Watersheds and 
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the shortest paths on critical percolation clusters with k = 1.04±0.02 |100) . are examples 
of unusual universality classes. Watersheds can also be dehned in higher dimensions m- 
By raising the water level through the landscape, different lakes are gradually forming 
and start merging to each other to form larger and larger lakes. It is likely to expect 
that at a certain water level, a percolation transition of the lakes would happen to form 
a giant lake so that it touches the borders and the water can indeed flow out of the 
landscape. This problem is closely related to the retention capacity of the landscapes 
addressed in [5SII571 [TT5] . Whether the percolation transition is critical or not, is strongly 
dependent on the spatial behavior of the correlations between the height variables |119j . 
For self-afhne surfaces with positive Hurst exponent H |120L 1121111221112311124j where 
the correlation behaves like ~ (1 — with the spatial distance r, there will not be 
a genuine percolation transition, while for a long-range correlated surface in which the 
correlation decays with the distance as ^ with H > 0, the percolation transition 

is critical [Mj and corresponding critical exponents change with H [12511126L112711128j . 
The fractal dimension of the watersheds [Util I129| as well as various geometric features 
are dependent on H |1301 11311 11321 I133L I134j . It has been verified numerically that 
the duality relation, as a characteristic property of conformally invariant fractals [135j . 
also holds for the perimeter of the largest cluster in the full range of Hurst exponents 

maiss]. 

For real landscapes, in addition, percolation theory provides an interesting descrip¬ 
tion for the global topography of Earth. It is found in [2T] that a percolation transition 
occurs on Earth’s topography in which the present mean sea level is automatically singled 
out as a critical level in the model. This finding elucidates the origins of the appearance 
of ubiquitous scaling relations observed in the various terrestrial features on Earth. This 
transition is shown to be accompanied by a continental aggregation which sheds light on 
the possibility of the important role played by water during the long-range topographic 
evolutions. The criticality of the current sea level also justifies the appearance of the scale 
(and conformal) invariant features on Earth, e.g., the fractal rocky coastlines |102j . with 
an intriguing coincidence of the dominant 4/3 fractal dimension in the critical model. 
The geometrical irregularity of the shorelines actually helps damping the sea waves and 
decreasing the average wave amplitude. A simple model is accordingly presented in |5D], 
which produces a stationary artificial shoreline related to the percolation geometry. A 
practical application of the discovery of the conformal invariance in the statistical prop¬ 
erties of the shorelines is that it allows one to analytically predict the highly intermittent 
spatial distribution of the flux of pollutant diffusing ashore |102) . 

This review paper is organized as follows. Section [^presents additional basic proper¬ 
ties of percolation model. Formulation of the fractal structure and critical properties of 
the model is given based on the scaling theory. Section [^describes different modifications 
and variants of percolation models which are mostly developed recently. In Section]^ we 
focus on the geometrical properties of the percolation model in two dimensions. The 
theory of SLE and its analytical consequences for the percolation model are briefly re¬ 
viewed. We also show that how one can map the criticality of a 3D Ising model onto a 
2D cross section of the original model. Section outlines the statistical properties and 
the modeling of the seashores and watersheds. We explain how percolation theory can 
describe the origins of the appearance of various fractal patterns on Earth. Finally, we 
briefly summarize our results in Section 
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2. Basic properties of the percolation model 


As we saw in the Introduction, percolation theory is concerned with the clustering 
properties of identical objects which are randomly and uniformly distributed through 
space with an occupation probability p. It is so simply defined yet so full of fascinating 
results. Although it is purely geometrical in nature, it embodies many of the important 
concepts of critical phenomena. In order to formulate the model, let us now define a 
bunch of geometric observables. For now we consider a system of infinite lattice size and 
then in Sect ion [2731 address the situation of finite size effects. For a given value of 

occupancy p, the nature of percolation is related to the properties of the occupied clusters. 
Two objects (e.g., occupied sites or bonds) belong to the same cluster if they are linked 
by a path of nearest-neighbor bonds joining them (this definition is slightly different 
for site and bond percolation). The connectedness is the essential characteristic of the 
percolation model denoted by the spanning probability Ps'. In the limit of an infinite 
lattice there exists a well-defined threshold probability Pc above which there suddenly 
emerges an infinitely large cluster that spans the system (Fig. |^. Therefore, for p > Pc 
there exists almost surely (i.e., Pg = 1) one infinite cluster of strength Poo{p) which 
denotes the probability for a given site to belong to the infinite cluster. All other clusters 
have finite size s at any arbitrary p, described by the cluster size distribution ns{p) i.e., 
the number of finite clusters of s connected sites, per lattice site. The probability that 
an arbitrary site belongs to a finite cluster of size s is then given by sns{p). Thus, the 
sum of all the probabilities that a given site belongs to either a finite size cluster or the 
infinite cluster must equal to p, i.e., 

'^sns{p) + P^{p)=p. ( 1 ) 


In the subcritical region, p < pc, Poo{p) is identically zero, and in the supercritical case 
p > Pc, Pcaip) is positive. Consequently, in the vicinity of the percolation threshold, the 
function Poo{p) is nonanalytic (Fig. [^. 

Now we need to define a quantity to measure the cluster size. Intuitively, at a low value 
of p cluster sizes are small and increase with p until the threshold where the spanning 
cluster dominates and is infinite in size, and thus the cluster size must diverge. Above 
the threshold we should remove the infinite cluster from our calculations otherwise it 
will always dominate. As the clusters get absorbed into the spanning cluster the typical 
size of those left goes back down again. Therefore, we should have a cluster size which 
increases, diverges at the threshold and then decreases again. 

By Bayes’ theorem, the probability a site belongs to a cluster of finite mass s given that 
the site is occupied is 


Ws{p) 


srigip) 

Y^gSUsipY 


( 2 ) 


Thus we can define the mean cluster size x{p) follows. 


X{p) = '^swg{p) 

S 


YgS^rig{p) 
YgSngip) ■ 


(3) 


This definition ([^ however, does not have information about the structure of the clusters 
e.g., their compactness and spatial extent. For each cluster of mass s one can instead 
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measure its radius of gyration Rg^s, defined by i?g ,, = (l/2s^) ^{ri — where the 
sum is over pairs of points on the cluster. For a given cluster of size s, the smaller radius 
of gyration is indicative of its higher compactness and lower spatial extent. 

The probability to find a finite cluster of large size s at a given point decreases expo¬ 
nentially with s in the subcritical regime m\- More precisely, there exists k{p) > 0, so 
that k(j>) —)■ oo as p —)■ 0 and k{p = Pc) = 0 such that 

Wsip) ~ s —oo. (4) 

It can also be shown that in the supercritical regime, the tail of the finite cluster size 
distribution has a rather smoother decaying form. In other words, there exist functions 
Ki(p) and K 2 (p), satisfying 0 < K 2 {p) < ki(p) < oo, such that [ld7) 

exp ( - < w,{p) < exp ( - (5) 

Note that the power is the order of the surface area of the sphere in d dimensions 

with volume s. This implies that the clusters are compact in the supercritical region |138j . 
At the critical point, this probability has a power-law decaying form of Ws{pc) ~ 
with some critical exponent d = 5{d) > 0—see also Fig. 



We can also define another length scale which is different from the average radius 
of clusters, i.e., the correlation length defined by the two-point correlation function 
gdf)- This is the probability that if one point is in a finite cluster then another point a 
distance r away is in the same cluster. This function then typically has an exponential 
decay given by a correlation length ^ 

9c{r) ~ r — oo. (6) 

The correlation length is a characteristic size of the cluster distribution which yields 
a maximum size above which the clusters are exponentially scarce. It is also the up¬ 
per bound of the scaling region where percolation clusters have a self-similar behavior. 
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Therefore we may define the correlation length ^ as an average distance of two points 
belonging to the same cluster 


.2 . 

Y.r9c{r) 


For a given cluster of mass s, one may replace in the above summation by the average 
squared distance between two cluster points i.e., 2i?g Moreover, with the probability 
sUs, a point belongs to an s-cluster, and since it is then connected to s sites, one may 
also replace gdr) by s^rig, giving rise to the following relation for the squared correlation 
length 1^ 


eip) = 






■ns{p) 




s{p) 


( 8 ) 


The above definitions of different characteristic observables are valid in all dimensions d. 


2.1. Percolation in d-dimensions 

2.1.1. Percolation on 

Let us consider the percolation problem on a hypercubic lattice 1/“ in d-dimensions. In 
d = I, it is a trivial task to find that the critical threshold should be Pc = 1 both for bond 
or site percolation models. For d > 1, there exists a critical threshold 0 < Pc < 1 below 
which all open clusters are finite and there is, almost surely, no infinite open cluster, 
and above Pc there exists an infinite open cluster with probability 1. It is rigorously 
known that no infinite open cluster exists at p = pc for d = 2 and d > 19 [52]. For other 
dimensions it is also conjectured to be held but its proof is viewed as an important open 
mathematical problem in the field. It is also known that for an infinite connected graph 
with maximum finite vertex degree A, the bond and site critical thresholds i.e., and 
p®, respectively, satisfy the following inequality |137j 

-^<pl<pl<l-{l-p\E- (9) 

In particular, p\ < p® where the strict inequality holds for a broad family of graphs. It 
can also be shown that for percolation on 1/ the percolation probability always satisfies 
P^^^\p) > Px\p) for allp and d, and consequently we have Pc(d-l-l) < Pc(d) |139lll40j . 

In order to further figure out the global feature of percolation, a natural question 
would be concerned about the number of possible infinite clusters which can coexist. It 
was shown by Newman and Schulman [141] that for periodic graphs at any arbitrary 
p, exactly one of the following three situations prevails with probability 1: The number 
of infinite open clusters can be either 0, 1 or oo. It has been proved in [142] that the 
third situation is impossible on l/. It has also been proved in [143] that there cannot be 


^The mean-square distance between two sites on an s-cluster is related to Rg.a, = 

2i?^ g, where the factor of 2 comes from counting each pair twice. The correlation length ^ can be defined 
as an average distance between two cluster sites. Whereas 2ijg_„ is the mean-square distance between 
two sites on an s-cluster, is this same distance averaged over all finite sizes s. The probability that a 
site belongs to s-cluster is sn^. There are s sites in each s-cluster. One can thus weight the average of 
2ijg g by s ■ sria, to obtain the squared correlation length Eq. Im. 
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Figure 7: Part of a Bethe lattice with coordination number 2 : = 3 IMI. The lattice sites are represented 
by open small circles at different shells k = 0^1,2, - . 


infinitely many infinite open clusters on any amenable graph Nevertheless there are 
some graphs, such as regular trees, on which infinitely many infinite clusters can coexist 

mi- 


2.1.2. Percolation on Bethe lattices 

Due to its distinctive topological structure, several statistical models even with in¬ 
teractions defined on the Bethe lattice |146) are exactly solvable and computationally 
inexpensive m- Various systems including magnetic models m and percolation 
|1481114H| , have been studied on the Bethe lattice whose analytic results gave important 
physical insights to subsequent developments of the corresponding research fields. The 
Bethe lattice is defined as a graph of infinite points each connected to z neighbors (the 
coordination number) such that no closed loops exist in the geometry—see Fig. A 
finite type of the graph with boundary is also known as a Cayley tree and possesses the 
features of both one and infinite dimensions: since Nk, the total number of sites in a 
Bethe lattice with k shells, is given as Nk = [z{z — 1)^ — 2]/(z — 2), the lattice dimension 


^For a finite set of vertices V whose edge boundary is denoted by OeV ={{u,v):u £ V, v ^ V}, 
the notion of amenability is related to whether the size of OeV is of equal order as that of V, or is 

much smaller. Denoting the Cheeqer constant of a graph G by h(G) = inf that is 

VCV(G):\V\<oo 1^1 

the minimal ratio of boundary to bulk of its nontrivial subgraphs, a graph is called amenable when 
h{G) = 0, and non-amenable otherwise. The simplest example of a non-amenable graph is the Bethe 
lattice with z > 3 for which the Cheeger constant is h = z — 2 | 140l I144| . It is also shown that 
Pc(G) < l/{h{G) + 1), so that for every non-amenable graph Pc{G) < 1. 
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defined by d = limfe_>,oo[lnfVfe/lnfc] is infinite. It is therefore often mentioned in the 
literature that the Bethe lattice describes the infinite-dimensional limit of a hypercubic 
lattice 1/ (It is also clear that a Bethe lattice with z = 2 is isomorphic to the positive 
ID lattice Z+). The Bethe lattice is indeed a very important substrate or medium on 
which the mean-field theories for various physical models become exact. 

As the number of shells grows the number of sites in the surface, or the last shell, grows 
exponentially z{z— 1)^“^. Therefore, as k tends to infinity, the proportion of surface sites 
tends to {z — 2)l{z — 1). By surface boundary we mean the set of sites of coordination 
number unity, the interior sites all have a coordination number z. Thus the vertices of a 
Bethe lattice can be grouped into shells as functions of the distances k from the central 
vertex. 

It is possible to show that the critical threshold for the Bethe lattice is Pc = l /(-2 — I) 
for any z > 3. Moreover, Poo = 0 indicates the subcritical region p < pc, while Poo > 0 
indicates the supercritical region p > pc in which Poo (p) is strictly increasing function of 
p. It is also known that at the critical point Poo(p = Pc) = 0, meaning that there is no 
infinite cluster almost surely at pc. 

Since the Bethe lattice has a tree structure, the only way to connect a point sitting at 
fcth shell to the origin is by a path of k edges. Thus the two-point correlation function is 
given by gdk) = p^, which decays exponentially fast as fc —>■ oo for all p < 1. Comparing 
this with relation it turns out that the correlation length should be ^(p) ~ — 1 /lnp. 
Therefore, unlike for the lattice, the correlation length on the Bethe lattice is finite 
for all 0 < p < 1 . 

The mean cluster size on the Bethe lattice can also be explicitly computed as follows. 


xip) = I 

( oo 


P < Pc, 
otherwise. 


( 10 ) 


In the next section we will see how x(p) —>■ oo as p approaches Pc from the left. 

A much wider class of interesting graphs is that of Cayley graphs (which includes 
also Cayley trees as a particular kind) of infinite, finitely generated groups. There, it has 
been shown m that the number JVoo of infinite clusters satisfies 

0 if pG[0,Pc), 

00 if P&iPcPu), (11) 

1 if pG(p„,l]. 

The parameter space is thus split into three qualitative intervals separated by two critical 
values Pc and p„. Some of intervals may be however degenerate or empty e.g., for Z"^ we 
have Pc = Pu, and for trees we have p^ = 1. 



2.1.3. Percolation on random graphs and networks 

The field of random graphs was started in 1959 by Erdos and Renyi |152Lll531ll54Lll55j 
whose work instigated a great amount of research in the field |15611157) . Random graphs 
have been extensively used as a probabilistic approach to study complex networks |158j . 
Many real-world complex networks such as the Internet |159j . social networks |160j . 
disease modeling |161j etc., share similar features e.g., they are large, sparse, scale-free, 
small worlds and highly clustered m- Consider a graph G„, with n number of vertices. 

(n) 

Denote the proportion of vertices with degree k in G„ by random variable . We 
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first call a random graph process sparse when lim„_>,oo , for some deterministic 

limiting probability distribution {pk}^o- Also, since {pk}^^o sums up to one, for large 
n, most of the vertices have a bounded degree, which explains the phrase sparse random 
graphs. We further call a random graph process scale-free with exponent r when it is 
sparse and when 


logPfc 

log{l/k) 


( 12 ) 


exists. This means that the number of vertices with degree k is proportional to 
an inverse power of k, i.e., Nk ~ where c„ is some normalizing constant. The 

requirement Nk = n then makes it reasonable to assume that r > 1 . 

Let the random variable denote the typical distance between two uniformly chosen 
connected vertices in G„. Then, we say that the random graph process is a small world 
when there exists a constant K such that with probability one, < K log n, in the limit 
of n —)■ cx). For ultra-small world random graphs we have almost surely < K log log n. 

The simplest imaginable random graph is the Erdos-Renyi random graph ER„(p), 
which arises by taking n vertices, and placing an edge independently between any pair 
of distinct vertices with some fixed probability p. This random graph is shown to exhibit 
a percolation phase transition in the size of the maximal component, as well as in the 
connectivity of the arising random graph. The phase transition in ER„(p), refers to a 
sharp transition in the largest connected component and serves as the mean-field case of 
percolation. 

The degree of a vertex in ER„ (j>) has a binomial distribution with parameters n and 
p = X/n. It is well known that for n is large, the proportion of vertices with degree k 
converges in probability to the Poisson distribution with parameter A i.e.. 


Pfc(A) 



as n 


oo, 


A: = 0,1,2,-. 


(13) 


In particular, po = is the fraction of isolated vertices. The control parameter A may 
thus be seen as the average degree of nodes. 

Therefore, ER„(p) is a sparse random graph process but not scale-free, since clearly it 
does not have a power-law degree sequence. However, in order to adapt the random graph 
to complex networks, we can make these degrees scale-free in a generalized random graph 
by taking the parameter A to be a random variable with a power-law distribution 
Erdos and Renyi have shown m that there is a critical value Ac = 1 below which 
ER„(p) random graph has almost surely no connected components of size larger than 
0{\ogn). At A = Ac, the graph has a largest component of size while for 

A > Ac, a drastic separation takes place between the largest cluster and all other smaller 
components and there appears a unique giant component of 0[n) containing a positive 
fraction of the vertices. 

The average fraction M[X) of clusters (number of clusters per vertex) can be shown |164j 
to be given by 

M(\) = -^(1 - PooiXf) + (1 - Poo(A))[l - In (1 - Poo(A))], (14) 

where Poo (A) denotes for the fraction of vertices in the largest component. 
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As mentioned earlier, the random graphs like ER„(p) are quite unlike real-world 
networks, which often possess power-law or other highly skewed degree distributions. The 
study of percolation on graphs with completely general degree distribution is, however, 
presented in [ 2 d) , giving exact solutions for a variety of cases, including site percolation, 
bond percolation, and models in which occupation probabilities depend on vertex degree. 

2.2. Percolation at and near criticality 

We have learned that the behavior of percolation process depends strongly on whether 
we are in the subcritical p < Pc or supercritical p > Pc regime. In the former case, all 
clusters are finite and their size distribution has a tail which decays exponentially—see 
Q. In the latter supercritical regime, there exists an infinite cluster with probability 
one and, the size distribution of other finite clusters has a tail which decays slower than 
exponentially—see ([^. At the vicinity of the critical point p ~ pc, however, there occurs 
an interesting phenomena characterized by a nonanalytic behavior of the order parameter 
Pooip) along with the divergent asymptotic behavior of the correlation length ^(p) and 
the mean cluster size xip) when p approaches Pc. Despite lacking their mathematically 
rigorous justifications, renormalization and scaling theory have made remarkable predic¬ 
tions about the behavior of the percolation problem near and at the critical threshold. 

2.2.1. Scaling hypotheses and upper critical dimension 

According to the scaling hypotheses it is possible to state, in general, the following 
relation for the number Ug (p) of hnite clusters per site 

ns(p) oc s~'^F(c[p)s'), s —oo (15) 

where r is a free exponent and E is a scaling function. Near the percolation threshold, c(p) 
is allowed to behave as a general power-law c(p) oc |p —Pc^^^, where a is another critical 
exponent. We can consider in general, the mth moment of the cluster size distribution 
defined by Mm{p) = '^aS'^ns{p) with m > 1. The following scaling relations are also 
conjectured to hold near the percolation threshold with different critical exponents j5^ 7 , 
a, A and ly 


percolation strength (p > pc): 

mean cluster size: 

mean cluster number per site: 

cluster moments ratio (m > 2 ): 

correlation length: 


Poo(p) ~P(p-Pc)^, 

(16a) 

x(p) ~r=^|p-pcr^. 

(16b) 

nc(p) ~A=^|P-Pcp"“, 

(16c) 

P=|-^ 

Mm{p) 

(16d) 

?(P) ^f'^\p-Pc\~'', 

(16e) 


which also define the critical amplitudes whose superscripts -I- or — refer to Pc being 
approached from above or below, respectively Universal combinations of these am¬ 
plitudes represent the canonical way of encoding the universal information about the 
approach to criticality [16511166] . While critical exponents can be determined working at 


®Eq. 1 16c I is valid only for the nonanalytic part of nc{p). 
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criticality, amplitude ratios also characterize the scaling region around the critical point 
which carry independent information about the universality class |167j . 

The probability for two sites separated by a distance r to belong to the same cluster also 
has a power-law decaying form gdr) — for large distances at p = Pc, introducing 

the anomalous dimension rj. 

These exponents, however, are not independent of each other but satisfy two sets of 
scaling and hyperscaling relations. The scaling relations can be easily read as 7 = v{2—rj), 
2 — a = 7 -I- 2/3 = (r — 1)/(T and /3 = A(t — 2). The hyperscaling relation, on the other 
hand, involves the number d of dimensions dv = 2 — a, and believed to be valid only 
for d < dc, where dc is called upper critical dimension. It is believed that when d > dc, 
the percolation process behaves roughly in the same manner as percolation on an infinite 
regular tree and their critical exponents take on the corresponding values given by mean- 
field theory: 

a =-1, 13 = 1, 7=1, = 5 = 2, A = 2, p = 0, ^=^- 

If these values are attained by percolation, then the hyperscaling relation gives dc = 6. 
The critical exponents are known exactly only in 2D, and not much is known rigorously 
in the general case. Critical exponents are universal in the sense that they depend only 
on dimensionality d, and not otherwise upon the individual structure of the underly¬ 
ing lattice. Theory of renormalization group (RG) lends support to the hypothesis of 
universality. 

2.2.2. Real-space renormalization group 

The Kadanoff picture of RG, called real-space RG, is based on coarse graining and 
rescaling procedure in which the lattice is iteratively divided into blocks of linear size b 
and then rescaled. When this scale transformation is iterated many times, RG leads to a 
certain number of fixed points. The fixed point equation for the rescaling transformation 
^ in percolation, has two solutions only: ^ = 0 , 00 . These are associated with 

the solutions to the fixed point equation in p-space, Tt,(p) = p, that is, p = 0,1 and Pc, 
representing the trivially self-similar states of the empty and fully occupied lattice and 
the nontrivial self-similar state at p = Pc, respectively. The critical exponent iz can then 
be given by 

^ ^ log(^) _ log(^) 

/ dTb{p) I \ / dRbjp ) , \ 

dp [ dp 

where the rescaling transformation Tb has been substituted with a real-space renormalisa¬ 
tion transformation Rb incorporating coarsening with rescaling. Note that Pc is identihed 
with p*, the nontrivial solution to the fixed point equation Rb{p*) = p*. The real-space 
renormalisation transformation Rb{p) is often chosen to be the probability of having a 
spanning cluster (or a majority of sites occupied) in the block. 

One can therefore summarize the real-space renormalisation transformation procedure as 
follows: i) Divide the lattice into blocks of linear size b. ii) Replace all sites in a block by 
a single block of size b occupied with probability Rb{p) according to the coarse graining 
procedure, iii) Rescale all length scales by the factor b. 

As critical exponents are determined by the large scale behavior, they are universal and 
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insensitive to details of lattice structure. However, real-space renormalization gives lim¬ 
ited results and does not give the critical exponents, except in the case of the Bethe 
lattice. 


2.3. Fractal structure of the critical percolation clusters 

Scaling theory asserts that whenever system is viewed on length scales smaller than 
the correlation length it behaves as it does at the threshold. At the critical point, ^ as 
the only length scale dominating the critical behavior of an infinite lattice, is diverging 
i.e., ^ —>■ oo. Disappearance of this scale at p = Pc, is reminiscent of scale invariance 
which implies the emergence of self-similarity in the geometric feature of the percolation 
clusters. The fractal properties persist even for p ^ pc with finite whenever the length 
scale to investigate the system is less than once it exceeds the geometry becomes 
Euclidean. Consider a percolation system which is viewed through a hypercubic window 
of size L‘^ where L ^ is the linear window size or can be regarded as the size of a 
finite system. Scale invariance then requires that the mean mass M of the cluster within 
the window would increase as a power-law with size i.e., M(^, L) ^ L'^f, where is the 
fractal dimension of the cluster. Above Pc on length scales L ^ the infinite cluster 
can be regarded as a homogeneous object which is composed of many cells of size 
i.e., M{^,L) ~ f {L/. These can be mathematically summarized in terms of the 
crossover function m as follows. 




where 


m(L/0 = 


constant for L <C 

for L > C. 


(17) 


The mass M, on the other hand, is proportional to L“Poo- Equating this with ( |17[ ) and 
rewriting in terms of (p—pc) using (16a I and (16e), yields the scaling relation d'j = d—(3/v. 
Since the exponents /? and v are universal, the fractal dimension df is also universal. 
The values of dj are known exactly only in 2D and d > dc = 6 as dj = 91/48 and 
4, respectively. In other dimensions the estimates exist only by numerical simulations. 
However the fractal dimension d'j by itself is not enough to characterize the percolation 
cluster. For example, a critical percolation cluster and a pattern of diffusion-limited 
aggregation (DLA) in 3D, share the same value of fractal dimension d^ ~ 2.5, while their 
fractal structures have completely different features. For its better characterization one 
can define the shortest path between two cluster points. The shortest path of length Imin 

jmin 

is shown to be self-similar which satisfies the scaling relation Imin ~ R ^ where R is 
the linear distance between the two points and, dj*™ denotes for the shortest path fractal 


dimension. The chemical dimension d'^^ is then defined by M ~ Imin ~ R‘^^ which 
means dj = d™*"dj^. The fractal dimension d™“ is known exactly only for dimensions 
d > 6 to be d™“ = 2. Even in 2D, despite its relevance, the fractal dimension of the 
shortest path is among the few critical exponents that are not known exactly |168L I169j . 
Now d™*" distinguishes between percolation cluster and DLA in 3D, by d™“ ~ 1.38 and 
1 , respectively. 

A fractal percolation cluster is composed of several other fractal substructures including 
its perimeter (hull), external perimeter, backbone, dangling ends and red sites (bonds), 
etc. For instance, the mean number of red bonds varies with p as A/,(p) ~ (p~Pc)~^i 
implying Nr which gives the fractal dimension of the red bonds d^ = 1/v, valid 
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in all dimensions [170| . 

In 2D, the perimeter and the external perimeter have the fractal dimension = 1 jA, and 
dj ^ = 4/3, respectively, which belong to the family of conformally invariant curves called 
SLE,^, with K = 6 and R = 8/3, respectively, satisfying the duality relation kR = 16. 

The fractal dimensions for the percolating cluster at criticality in ER„(p) random 
graph and random scale-free networks with n equal to the number of vertices and degree 
distribution pk = Cnk~'^, are also reported in m- There has been shown that the 
fractal dimension of the spanning cluster is dj = 4 for r > 4 and d^ = 2(r — 2)/(t — 3) 
for 3 < T < 4. Note that the result for r > 4 is in agreement with that of the regular 
infinite dimensional percolation. As we mentioned earlier in subsection |2.1.3[ on a random 
network in the well connected regime, the average distance between sites is of the order 
n, and becomes even smaller in ultra-small world random graphs. However as 
discussed above, for d > 6 we have dj = 4 and d™*" = 2, consequently d°^ = 2. 
Therefore the average chemical distance Imin between pairs of sites on the spanning 
cluster for random graphs and scale-free networks with r > 4 at criticality behaves as 
Imin ~ VM, i.e., the distances become much larger at criticality. 


3. Variants of percolation 

In the previous section, we studied the standard version of percolation i.e, Bernoulli 
percolation, in which all random occupations of bonds or sites take place irreversibly and 
independently on either a Euclidean lattice or a random graph. Over the past decades 
and years, numerous variations and modifications of the percolation model with a huge 
variety of applications in many fields were introduced. An important development first 
established by Fortuin and Kasteleyn, is the connection between bond percolation and a 
lattice statistical model. This consideration has led to a formulation of the percolation 
problem as a limiting case of the general Potts model, which was extremely useful, for 
many of the techniques readily available in statistical mechanics has been applied to per¬ 
colation m- Many variations were motivated by consideration of spatial correlations, 
anisotropy, nonlocality, explosivity etc. ’Explosivity’ in particular, deals with the search 
for those perturbations of cluster-merging rules which change the order of percolation 
transition from the second order (continuous) to the first order (discontinuous). This is 
known as explosive percolation in which a macroscopic connected component emerges in 
a number of steps that is much smaller than the system size. This has recently become 
a subject of enormous interest m, including openings toward other subjects such as 
jamming in the Internet m, synchronization phenomena |I75L 1176) and analysis of 
real-world networks Let us first start by introducing the family of explosive perco¬ 

lation models and then turn our attention to some other variants and modifications. The 
number of papers appeared on explosive percolation was itself explosive and we cannot 
cover all subjects and ideas of course, we rather try to outline some original and selective 
topics. 


3.1. Explosive percolation 

3.1.1. Achlioptas process 
Having introduced in subsection 

briefly ER„ model), is composed of n isolated vertices whose each pair of vertices is 


2.1.3 the classic Erdds-Renyi random graph (or 
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chosen uniformly at random in each step, and connected by an edge {ei}. At any given 
moment, a cluster is defined as the set of vertices each of which can be reached from any other 
vertex in it by traversing edges. If tn denotes for the number of added edges at time t, it is 
known that the fraction of vertices in the largest cluster undergoes a continuous phase transition 
at the critical time tc = 1/2. 

At a Fields Institute workshop in 2000, Dimitris Achlioptas suggested a class of variants 
of the classical ER„ model where a nonrandom selection rule, the so-called Achlioptas process, 
is additionally imposed which tends to the delay (or acceleration) in the formation of a large 
percolating cluster. This has then received much attention in recent years |178l 1179111801 [181] . 
Concretely, consider a model that, like ER„, starts with n isolated vertices and add edges one 
by one. The difference is that to add a single edge, first two random edges {ei,e 2 } are chosen, 
rather than one, each edge is chosen exactly as in ER„ and independently of the other. Of these, 
only one should be selected according to the selection rule, and then inserted in the graph. The 
other edge is discarded. Clearly, if one always resorts to randomness for selecting between the 
two edges, the original ER„ model is recovered. 

Achlioptas originally asked if it is possible to shift the critical point of this phase transition 
by following an appropriate selection rule. One rule that can naturally be imagined is the 
product rule: Of the given potential edges, pick the one which minimizes the product of the 
sizes of the components containing the four end points of { 01 , 02 }. This rule was suggested in 
|182| as the most likely to delay the critical point. Another rule is the sum rule, where the size 
of the new component formed is minimized. 

A selection rule can be classified as a bounded-size or an unbounded-size rule. In a bounded- 
size selection rule, decisions depend only on the sizes of the components and, moreover, all sizes 
greater than some (rule-specific) constant K are treated identically. For example A = 1 is the 
Bohman-Frieze (BF) rule, where ei is chosen if it joins two isolated vertices, and 62 otherwise. 
It was rigorously proven in | 178 | for a much simpler rule, that such rules are capable to shift 
the threshold. Moreover, the percolation transition is strongly conjectured to be continuous for 
all bounded-size rules For unbounded-size rules in contrast, extensive simulations |58| 

strongly suggested that, the product rule in particular, shows much more interesting behavior 
than just shifting the critical point, it exhibits a discontinuous percolation transition. The 
numerical evidence showed that the fraction of vertices in the largest cluster jumps from being 
a vanishing fraction of all vertices to a majority of them instantaneously, i.e., the largest cluster 
grows from size at most y/n to size at least n/2 in at most steps, that is, at the phase 

transition a constant fraction of the vertices is accumulated into a single giant cluster within a 
sublinear number of steps. 

However, this conjecture has been rigorously disapproved by Riordan and Warnke [7211183] 
in 2011, by showing that it cannot be a discontinuous transition, but a continuous one. In fact, 
their argument shows continuous phase transitions for an even larger class of processes |183| 
called Z—vertex rules (every Achlioptas process is a 4—vertex rule). Their results state that 
continuity of the phase transition is such a robust feature of the basic model that it survives 
under a wide range of deformations and thus all Achlioptas processes have a continuous phase 
transition. This however neither means that the connectivity transition at pc is characterized 
by a usual power-law divergence of the order parameter |184l 118511186 ] nor that it cannot be 
followed by multiple discontinuous transitions. In particular, Nagler et al. have proven that 
for certain Z—vertex rules the continuous phase transition can have the shape of an incomplete 
devil’s staircase with discontinuous steps in arbitrary vicinity of Pc [Ml- 

Therefore, heuristics or extrapolations from simulations suggesting explosive percolation in 
mean-field models seem to be wrong in the scaling limit. Nevertheless, a discontinuous transition 
is possible when one departs far enough from the ER„ model. 
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Figure 8: Illustration of the first four steps of time evolution of a graph Gn with n = 8 number of 
vertices, according to the half-restricted process. The vertices surrounded by the dashed line belong to 
the restricted vertex set Rf(G) with / = 0.5. 


3.1.2. Half-restricted process 

Half-restricted process |188| is a variant of the Erdos-Renyi process which exhibits a discon¬ 
tinuous phase transition. Intuitively, a discontinuous transition can only occur if the growth of 
large clusters is suppressed and clusters with medium sizes become abundant. After a number 
of steps, such medium-sized clusters merge and a giant cluster emerges drastically. As we saw 
in the last subsection |3.1.1[ the idea to do this in an Achlioptas process was to select an edge 
that connects smaller components, which did not lead to an explosive percolation. A different 
approach is however considered in the half-restricted process: In each step, two vertices are 
connected by an edge, but one of them is restricted to be within the smaller components. 

Consider a graph Gn including n labeled vertices vi,V 2 ,... ,Vn sorted ascending in the size of 
the clusters that they reside in. Vertices with the same cluster size are sorted lexicographically. 
For a given value of 0 < / < 1, a restricted vertex set Rf{G) is defined which is composed 
of the [/nj number of vertices in smallest components {\_fn\ means the floor of fn, i.e., the 
largest integer less than or equal to fn). At each time step t > 1, one edge is added between 
two chosen vertices Vr and Vu. Starting with an empty graph at t = 1, at every step t, first 
the restricted vertex set Rf{G) is recognized, and one restricted vertex Vr, which belongs to 
Rf{G), is uniformly chosen at random. Then independently, one unrestricted vertex Vu is 
chosen uniformly at random from the whole vertex set Gn. An edge is then added between 
two vertices Vr and Vu if is not already present. Figure [^illustrates the first four steps of time 
evolution of an Gn graph with n = 8 and / = 0.5, according to the half-restricted process. 

It has been shown in | 188 | that, unlike the Achlioptas process, the half-restricted process 
exhibits explosive percolation with a discontinuous phase transition for any value of / < 1. 
Although the evolution of the largest cluster over the first n edges in the half-restricted process 
may behave very similar to that of the Achlioptas process | 188| , they are fundamentally different 
in the nature of phase transition as well as mathematical structure. 

3.1.3. Spanning cluster-avoiding process 

For the continuous phase transitions, the Erdos-Renyi process on random graphs accounts 
for the mean-field limit of the standard percolation on a Euclidean lattice. Explosive percolation 
models in Euclidean space have also been extensively studied |6UI I73 | , and the numerical results 
suggest discontinuous transitions |189| . However, due to the lack of analytic arguments, the 
order of explosive transition in Euclidean space is still elusive. This would therefore be of 
interest to clarify the order of the explosive transition in Euclidean space and on random graphs 
in a unified manner. To this aim, a model called the spanning cluster-avoiding (SCA) model was 
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introduced [nsmss]. The spanning cluster in Euclidean space actually plays the same role as the 
giant cluster (or component) in random graph models. The SCA model starts by considering a 
finite hypercubic lattice Z'* in d dimensions of size N and unoccupied bonds. Inspired by the 
best-of-m model at each time step t, number of m unoccupied bonds are chosen randomly 

and classified into two types: bridge and nonbridge bonds. Bridge bonds are those that upon 
occupation a spanning cluster is formed. SCA model avoids bridge bonds to be occupied, and 
thus one of the nonbridge bonds is randomly selected and occupied. If the m potential bonds are 
all bridge bonds, then one of them is randomly chosen and occupied. Once a spanning cluster is 
formed, restrictions are no longer imposed on the occupation of bonds. This procedure continues 
until all bonds are occupied at t = 1. 

Extensive numerical simulations and theoretical results |190| have shown that the explosive 
transition in SCA model in the thermodynamic limit, can be either discontinuous or continuous 
for d < dc = 6 depending on the number of potential bonds m. In other words, there is a 
tricritical value mc(d) = d/(d — dsn) for d > dss, where dsB denotes for the fractal dimension 
of the set of bridge bonds [na. such that if m < me, the transition is continuous at a finite 
threshold and discontinuous, in the thermodynamic limit, for m > at the trivial percola¬ 
tion threshold at t = 1, when all bonds of the system are occupied. The formula for mc(d) is 
valid only for d < dc = 6 at which mc{dc = 6) = oo. 

For d > dc, i.e., in the mean-field limit, the transition is shown to be continuous for any finite 
and fixed value of m. However if m varies with the system size N, a discontinuous transition 
can also take place for d > dc. More precisely, there exists a characteristic value rric ~ In such 
that when m increases with N slower than rric, the transition is continuous in a finite critical 
time tc, and when m increases with N faster than rric, the transition is discontinuous at the 
trivial percolation threshold at f = 1. If m increases with N as rric ~ In N, then a discontinuous 
transition occurs at finite tem, which is neither nor unity. Some necessary conditions for a 
non-trivial and a trivial discontinuous percolation transition have been recently proposed | 193| . 

The idea to obtain a discontinuous percolation transition by controlling the largest cluster 
alone was actually addressed before m, in which an acceptance method was used to system¬ 
atically suppress the formation of a largest cluster on the lattice. There has been also shown 
that the cluster perimeters are fractal at the percolation threshold with a fractal dimension of 
1.23 ± 0.03, statistically indistinguishable from that of watersheds. 

3.2. Non-self-averaging percolation 

Fractional percolation, as a non-self-averaging percolation, was introduced | 194| to describe 
some of the main features of crackling noise as a characteristic feature of many systems when 
pushed slowly, e.g., the crumpling of paper m\, earthquakes and the magnetization 

of slowly magnetized magnets. In the fractional percolation the relative size of the largest 
component Smax, as the order parameter, exhibits many randomly distributed jumps after a 
critical threshold Pc whose discontinuities survive even in the thermodynamic limit. A fractional 
growth rule is used in this model which induces a certain type of size homophily among clusters: 
Connection of two clusters with a similar size is in priority, while the size of the larger cluster 
has already been rescaled by a target fraction factor 0 < / < 1. The model starts by considering 
a network of n isolated vertices with no edges. At each step, three different vertices Vi,V 2 ,V 3 are 
chosen uniformly at random, no matter if they are in the same cluster. Let Si > S 2 > S 3 denote 
the sizes of the (not necessarily distinct) clusters that they are contained in. Then an edge is 
added between the two vertices Vi and Vj for which Aij := fSi — 1 < i < j < 3 is minimal. 

If there are multiple options, one is randomly chosen. Moreover when only two clusters are left 
in the system they will be connected—see Fig. for further illustration. 

It has been shown that after the first transition, for p > pa and n —>■ 00 , the size of the largest 
component either stays constant or increases fractionally by at least a factor of r = //(I -|- /). 
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Figure 9: Illustration of the cluster merging process according to the fractional percolation with target 
fraction / = 1/2. The link between vertices vi and V 2 is established because A 12 = —5 is minimal. 


As the first transition is point-continuous [ziiinzi, the process necessarily undergoes infinitely 
many discontinuous transitions arbitrarily close to the first dynamical transition point pc, and 
the order parameter increases stepwise. The height St of the fcth step down the staircase has 
the form <5^ ~ (1 + giving rise to the scaling distribution of the jump size as ds ~ 

Both the size of the jumps, and the transition points are stochastic even in the thermodynamic 
limit. 

The process is therefore shown to be non-self-averaging in the sense that the relative variance 
of the size of the largest component given by {s^ax)/— 1, does not vanish for p > Pc 
as n —>■ 00 . This is in contrast to classical models (e.g., the ER„ model) for which the order 
parameter converges to a nonrandom function in the thermodynamic limit. 

The characteristics of fractional percolation are robust against an arbitrary variation of the 
target fraction factor /. Assuming a time-dependence form of f{k) = a/k where 0 < a < 1, 
one finds that the jump size distribution function decays faster than ~ characterized by the 
power-law fluctuations ds ~ s~^ with /3 = (1 + a)/a. 

For other models that exhibit non-self-averaging in percolation we refer to the literature 

[TggllTMlIMJll^ . 

3.3. Correlated percolation 

In the percolation models studied so far, all random occupations have been considered to be 
independent of each other with no spatial correlations. But this is not always the case when, for 
instance, percolation theory is applied to study transport and geometric properties of disordered 
systems [202112031120411205112061 120711208 | , since the presence of disorder usually introduces 
spatial correlations in the model. For sufficiently short-range correlations, the properties of 
the model will be the same as those of uncorrelated percolation. By increasing the range of 
correlations however, they may have a relevant contribution leading to new fixed points in 
the renormalization group study of the model. These can be quantified in terms of how the 
correlation function g(r) falls off at large distances r: When the correlations are short-range 
with a fall-off faster than r~'^, then according to the Harris criterion |209 |. they are relevant if 
di/ — 2 < 0, where n is the correlation length exponent for the uncorrelated percolation model. 
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Since for the percolation model we always have the hyperscaling relation dv — 2 = —a > 0, so 
short-range correlations do not change the critical behavior m\- For long-range corralations of 
the power-law form g{r) ~ with 2H < d instead, the extended Harris criterion | 210 | states 

that the correlations are relevant if Hv—1 < 0. In this case, the new correlation length exponent 
is given by the scaling relation i/h ~ 1/H. This relation has been then verified numerically in 
|M [20Hl[mim^[2T3| . Thus the critical exponents in a long-range correlated percolation can 
change depending on how the correlations decay with the spatial distance. 

Such a power-law decay of the spatial correlations is a typical characteristic feature of the 
height profile {h(x)} in random grown surfaces, where /i(x) is the height at the lattice site at 
position X. This indeed provides a convenient way to tackle with the correlated percolation on 
a lattice. For self-affine surfaces for which g{r) ~ (1 — r^^), it is shown in |119| that even in the 
thermodynamic limit, the percolation transition is only critical for H = 0. Percolation on these 
surfaces is actually governed by the largest wavelength of the height distribution, and thus the 
self-averaging breaks down. For long-range correlated surfaces where g{r) ~ r~^^, in contrast, 
the transition is critical and the self-averaging is recovered. Depending on the value of H, the 
correlation-length exponent is given as follows 


r 1/H if 0<H<l/u, 
\ y if H > l/v, 


(18) 


where v = 4/3 for a Euclidean lattice in 2D. It is thus natural to expect that other critical 
exponents also depend on the value of H. Such a dependence is numerically verified for the 
fractal dimensions of the largest cluster its perimeter d^jj and external perimeter df//, 

shortest path, backbone, and red sites [nsi. It has also been shown that, within the numerical 
accuracy, the hyperscaling relation d — (7/r + 2j3H)/r’H = Ih/vh + 2(d — dfu) is fulfilled by 
the exponents. The duality relation is numerically shown to be valid {df// — 
though the theoretical verification of these observations is still lacking. 


3.4- Bootstrap percolation 

The bootstrap percolation problem [214[I215] and its obvious variants deal with the dynamics 
of a system composed of highly coupled elements, each of which has a state that depends on 
those of its close neighbors. It has played a canonical role in description of a growing list of 
complex phenomena including crack propagation [216] . neuronal activity [ 2 T 7 I [2181 [219] . and 
magnetic systems | 220 | among others. 

The standard bootstrap percolation process on a lattice as the spread of activation or infection 
is defined according to the following rule with a given fixed parameter k > 2: Initially, each 
of the sites is randomly infected (or activated) with probability p and uninfected (inactivated) 
with probability 1—p, independently of the state of the other sites. Every infected site remains 
infected forever, while each uninfected one which has at least k infected neighbors becomes 
infected and remains so forever. This procedure is continued until the system reaches the stable 
configuration which does not change anymore i.e., when no uninfected site has k or more infected 
neighbors. The main question which arises is concerning the percolation of the infected cluster 
i.e., if there emerges a giant spanning cluster of infected sites of size 0{n) by the end of the 
process. 

Bootstrap percolation has been extensively studied on 2D and 3D lattices [22II1222112231 [224] 
(and references therein), including the proof of the existence of a sharp metastability threshold 
in 2D |221 | which has then been generalized to arbitrary d dimensional lattices |222[ 1223 ] . In 
particular, Schonmann [225] proved that on the infinite lattice Z'*, the percolation threshold 
Pc(Z‘^,fc) = 0 if fc < d, and pc(Z‘^,A:) = 1 otherwise. The finite size behavior (also known as 
metastability) was studied in [22411226[ I227| . and the threshold function was determined up to 
a constant factor, for all 2 < fc < d, by Cerf and Manzo [227]. The first sharp threshold was 
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determined by Holroyd mu, for fc = 2 on a finite 2D lattice of linear size L, who proved 
that 

TT^ 

Pc(Z, d = 2, A: = 2) ~ -, as L —>■ oo. (19) 

^ ^ ^ 181nL’ ^ ^ 

This has been recently |228| generalized to the finite lattice of linear size L as follows 

Pc{L, d,k) ~ ^ , for 2 < k < d, as L —>■ oo, (20) 

Viii(fc-i) ^J 


where In(fe) denotes for an fc—times iterated logarithm, Injj.) L = In (ln(fc_i) L). Althongh the 
fnnction X{d, k) is not exactly known, but it can be shown that it always has a finite value with 
the following properties |228l : A(2,2) = 7r^/18, A(d, 2) ~ (d — l)/2, and A(d, d) ~ rr^/Od as 
d —>■ 00. Clearly the result (191 is a special case of ( |20[ ) for d — 2 and k = 2. Surprisingly, 
there exist some predictions for the asymptotic threshold in 2D based on numerical simulations 
which differ greatly from the rigorous result ( |19[ ). For example, based on simulations in |229 |. 
the estimate Pc(L, d = 2, fc = 2) In L = 0.245 ± 0.015 is reported whereas the rigorous prediction 
is 7r^/18 = 0.548311 ■ • • . This apparent discrepancy between theory and experiment has been 
rigorously addressed in | 230 | to be due to a very slow convergence in the asymptotic limit when 
L —^ oo. 

Bootstrap percolation has also been studied on the random regular graph [23111232] , and 
infinite trees |233 | as well. In the context of real-world networks and in particular in social 
networks, a bootstrap percolation process can be imagined as a toy model for the spread of ideas 
or new trends within a set of individuals which form a network. In this sense, the bootstrap 
percolation has been recently | 234 | studied on the power-law random graphs of n vertices in 
which is shown that a giant cluster of size 0(n) emerges, with high probability, above a certain 
critical sublinear number of the initially infected nodes. This behavior is in sharp contrast 
with that observed [133 [13^ on an Erdos-Renyi random graph in which there will be no 
evolution when the number of initially infected vertices is sublinear. 

The final remark is regarding the close relation between the bootstrap percolation and 
another well-known problem in graph theory, that of the k-core of random graphs [37I[^E35| . 
The fc-core of a graph is the maximal subgraph for which all nodes have at least k neighbors 
within the fc-core. However there is a difference between the stationary state of the bootstrap 
percolation and the fc-core | 236| . Bootstrap percolation is an infection process, which starts from 
a subset of source nodes and spreads over a network according to the infection rules described 
earlier. The fc-core of the network can be found as an asymptotic structure obtained by a 
subsequent pruning of nodes which have less than fc neighbors. 


3.5. Directed percolation 

Considering a porous rock as a random medium in which neighboring pores are connected 
by small channels of varying permeability, an important problem in geology would be how 
deep the water can penetrate into it. Clearly the ordinary percolation model is not applicable 
to describe such phenomena since the gravity has weighted a specific direction in space i.e., 
the water propagation is not isotropic but directed. Directed percolation, introduced in m , is 
an anisotropic variant of standard isotropic percolation which introduces a specific direction 
in space. It accounts for one of the most prominent universality classes of nonequilibrium 
phase transitions, playing a similar role as the Ising universality class in equilibrium statistical 
mechanics |239 |. Directed percolation models exhibit a continuous phase transition with a 
fascinating property of robustness with respect to the microscopic dynamic rules. It turned 
out to describe a wide range of spreading models e.g., contact process [2401 1241] . epidemic 
spreading without immunization |242| , and forest fire models | 243[ 124411245] . For this model, 
the percolation can only occur along a given spatial direction. Regarding this direction as a 
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Figure 10: Directed bond percolation on a tilted square lattice started from a source active vertex Vs 
at t = 0. Each bond (and corresponding destination site) is activated with probability p, shown by 
solid lines with arrows (filled circles). The cluster of active sites connected by a directed path to Vs is 
indicated in blue (dark). The vertical direction corresponds to time, and the dashed lines identify the 
sets of available vertices at time t. 


temporal degree of freedom, directed percolation can then be viewed as a dynamical process 
in d + 1 dimensions. In this sense percolation on a dynamic network can be mapped onto the 
problem of directed percolation in infinite dimensions | 246 |. 

Starting from a source active (occupied or wet) vertex at t = 0 on a tilted hypercubic 
Z'* lattice, directed bond percolation, as a dynamical process, can be interpreted as follows. As 
illustrated in Fig. |10[ at the next time step, each of the downward bonds emanating from Vs is 
randomly occupied by an arrow with probability p which corresponds to the destination site to 
become active. This procedure is continued row by row until the system reaches an absorbing 
state, i.e., a configuration that the model can reach but it cannot escape from there. There 
exists a critical threshold pc, that for p < pc the average number of active sites grows for a short 
time and then decays exponentially. For p > pc there is a finite probability that the number of 
active sites diverges as t ^ oo. In this case activity spreads within a so-called spreading cone. 
At p = Pc which separates a non-fluctuating absorbing state from a fluctuating active phase, a 
critical cluster is generated from a single source whose scaling properties are characterized by 
a number of critical exponents. Note that, in directed percolation each source vertex generates 
an individual cluster, so the lattice in this case cannot be decomposed into disjoint clusters. 
Therefore, depending on the value of p, activity may either spread over the entire system or 
die out after some time. The latter case i.e., the absorbing state is a completely inactive state 
which can only be reached by system but not be left. Thus detailed balance is no longer obeyed 
and that is why the process is nonequilibrium. 

Despite of its very simple rules and robustness, its critical behavior is highly nontrivial and exact 
computation of the critical exponents seems impossible. Even in 1 -|- 1 dimensions, no analytical 
solution is known, suggesting that is a non-integrable process. This may be related to the fact 
that directed percolation, unlike the ordinary percolation, due to the lack of symmetry between 
space and time, is not conformally invariant. 

Although it is not an equilibrium model, the critical behavior of the directed percolation 
shares a very similar picture as in the ordinary percolation. A phenomenological scaling the- 


27 

















ory can be applied to describe its criticality. Considering the density p{t) of active sites as 
an order parameter of a spreading process, observations justify that in the active phase p{t) 
decays and eventually saturates at some stationary value Ps{t). Near the critical point, the 
stationary density is then turned out to satisfy a power-law relation p^ ~ {p — pc)^, where /3 is a 
universal critical exponent that only depends on dimensionality. The other important quantity 
is the correlation length whose definition needs a special care, since in this case, time is an 
additional dimension which should be distinguished from the spatial dimensions. Let us denote 
the temporal and spatial correlation lengths by the indices and respectively, which are 
independent of each other. Then close to the transition, these length scales are expected to 
diverge as ~ \p — Pc\‘'^^ and ~ — with generally different critical exponents and 

The two correlation lengths are related in the scaling regime by ~ ^x where 2 : = v\\/v_l 
denotes for the so-called dynamic exponent. In many models, the universality class is given by 
the three exponents /3, i/y and 

For d > dc = 4, the values of the critical exponents are believed to be given by the mean-field 
theory as /3 = 1, z^y = 1 and = 1/2 [239]. For d < dc however, there are no exact results 
neither for critical exponents nor thresholds. However, a very precise estimates of critical expo¬ 
nents and thresholds on several lattices can be found in [24711248 | in 1 -|-1 dimensions, and in a 
recent work [249] (and references therein) in higher dimensions. 

Let us conclude this subsection by stating the following conjecture first made by Janssen |25U| 
and Grassberger |251| inspired by the variety and robustness of directed percolation models. The 
statement is that any model which fulfills the following conditions should belong to the directed 
percolation universality class: i) The model exhibits a continuous phase transition into a unique 
absorbing state, ii) The transition is characterized by a positive one-component order parameter, 
iii) The dynamic rules involve only short-range interactions, and iv) The system has no special 
attributes e.g., additional symmetries or quenched randomness. Although this conjecture is not 
yet rigorously proven, it is highly supported by numerical evidence. 

4. Percolation in two dimensions 

At the critical point, the correlation length diverges and the system becomes scale invariant 
and thus scaling hypothesis applies. However the scaling hypothesis alone cannot determine 
the critical exponents though it may give some relations among them. Therefore one may call 
for some possible stronger symmetries e.g., conformal invariance, to determine the exponents. 
Of course scale invariance does not imply conformal invariance at least at the level of the 
superficial mathematical definition. How can thus conformal invariance to be enhanced from 
scale invariance? In two dimensions, one can rigorously show that scale invariance is enhanced 
to conformal invariance under the following assumptions [TnillHimj: i) unitarity, ii) Poincare 
invariance (causality), iii) discrete spectrum in scaling dimension, iv) existence of scale current, 
and v) unbroken scale invariance. Fortunately most interesting classes of 1 + 1 quantum field 
theories, as the scaling limit of various 2D lattice models in statistical mechanics, satisfy these 
assumptions. In two dimensions however the conformal symmetry is extremely powerful since 
conformal transformations correspond to analytic functions to be used in statistical mechanics 
to characterize universality classes. Under conformal transformations the lengths are rescaled 
non-uniformly while the angles between vectors are left unchanged. 

Many exact results have then been obtained for percolation model in two dimensions using 
methods of conformal field theory (CFT). Among them, Cardy’s conjectured formula [25411255] 
for the crossing probability is one of the famous ones. For a percolation model defined in a unit 
disc I 2 I < 1, the probability that there exists at least one cluster which contains at least one 
point from each of two disjoint intervals on the boundary of the disc whose ends are assigned by 
{zi,Z 2 ) and ( 21 , 22 ) respectively, has the following explicit form as a function of the cross-ratio 
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where 2F1 is the hypergeometric function. Ps{r]) is invariant under transformations of the unit 
disc which are conformal in its interior (but not necessarily on its boundary). Of course, this 
probability is interesting only at the critical point. For p < pc, since all clusters are finite, in the 
scaling limit we will have Ps = 0 , while for p > pc the infinite cluster always spans, so the limit 
is 1. This result is valid as long as there are only short-range correlations in the probability 
measure, independent of whether the microscopic model is formulated as bond, site, or any other 
type of percolation. 

Moreover it has been shown that the whole probability distribution of the total number of such 
distinct crossing clusters is conformally invariant in the scaling limit, depending only on 77 . In 
particular the mean number Afc of such crossing clusters is | 256 | 
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The other exact result obtained by using conformal field theory and Coulomb gas methods, 
is for the number Af{A) of percolation clusters of enclosed area greater than or equal to A at the 
critical point in two dimensions [257]. It has been shown to behave as the following power-law 
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with a proportionality constant C — l/iy/i-K that is universal. 

A rigorous proof of Cardy’s formula for site percolation on the triangular lattice was discov¬ 
ered m by Smirnov using advantages of the theory of stochastic Lowner evolution (SLE). 


4-1. Stochastic Lowner evolution 

In probability theory, the stochastic Lowner evolution (SLE^) with parameter k, introduced 
by Schramm |82| in 2000, is a family of random planar curves that have been proven to be 
the scaling limit of interfaces in a variety of 2D critical lattice models in statistical mechanics. 
Here we give a brief description of the so-called chordal SLE in its standard setup, and refer to 
(531 IMl 1251] for more detailed information about it, and other variants of SLE. 

Consider a curve 7 ( 1 ) that emanates at t = 0 from origin on the boundary of the upper 
half-plane H, and goes to infinity as f —>■ 00 . The curve might intersect but should not cross 
itself during the evolution. Let us define the hull Kt as the union of the curve and the set of 
points got trapped by the curve up to time t, which are not reachable from infinity without 
crossing the curve. According to the Riemann mapping theorem, there should be an analytic 
function (?t(«) which maps H\A't into the H itself (HI\Art is a simply connected domain that 
contains all those points of HI that are not in Kt). The map gtiz) can be uniquely determined 
by imposing the following hydrodynamic normalization at infinity, 

2t 

gt{z) = z + — -, as a —>■ 00 , (24) 

where the coefficient 2t is due to a conventional parameterization of 7 . It can be then shown 
that the evolution of the tip of the curve, can be given by the Lowner differential equation 

2 

gt{z) - Ct’ 

29 


dtgt{z) 


go{z) = z, 2 ; e H, 


(25) 










where Q is a continuous function but not necessarily differentiable. In order to have conformally 
invariant random curves which behave geometrically as they should to encode the statistics of 
critical interfaces, Schramm argues that |82| they should have two properties: Markov property, 
and stationarity of increments property. With these two properties and reflection symmetry, Q 
can only be proportional to a ID standard Brownian motion i.e., so that {C,t) = 0 and 

((Ct — Ct')^) = — The diffusivity k is the only parameter whose different values correspond 

to different universality classes of critical behavior. For example, k = 6 corresponds to the 
percolation universality class in which cluster boundaries in the continuum limit are described 
by SLEg. In fact SLE^ is in general conformally covariant under domain changes, and only for 
the special case of percolation with k = 6 is conformally invariant. 

For ft = 0, SLE curve is a vertical straight line, and when ft increases the curve become more 
rough with fractal dimension df = 1 + ft /8 for ft < 8 and 2 for ft > 8 [259]. For 0 < ft < 4, the 
SLE curve does not intersect itself or the real axis {dilute phase), while for ft > 4, it intersects 
but does not cross itself and the real axis on all length scales (dense phase). SLE is a space-filling 
curve when ft > 8 . In the dense phase with ft > 4, there is a duality conjecture stating that the 
exterior frontier of an SLE^ hull looks locally as SLE^ with ft = 16/ft |135l I26(J| . For example, 
the external perimeter of a percolation cluster in the continuum limit is believed to be described 
by SLEg /3 which is also expected to describe the scaling limit of planar self-avoiding random 
walk (SAW), although there is no complete mathematical proof yet |261| — see also Fig. |11| 



Figure 11: For a putative percolation cluster, the perimeter (or hull) is the union of the solid and 
dashed boundary lines which contains many fjords, and the external perimeter is the union of the solid 
and dotted boundary lines with fractal dimensions and , respectively. The external perimeter is 
obtained by closing off all narrow passageways. The duality relation states that (d^^ — l)(d^ ~ 1) = 1/4. 

In addition, SLE^ has two more special properties for k = 6 (locality) and k = 8/3 (restric¬ 
tion). Let D be a simply connected region in H connected to the real axis which is at some finite 
distance from the origin. Consider two SLE^ processes from the origin to infinity, one in the 
domain H and another in the domain BI\D. If these two processes have the same distribution 
up to the hitting time of the set D, then the SLE„ has the locality property. Such a property is 
expected for the percolation cluster boundaries with ft = 6 (Mum], and for no other values of 
ft. Moreover, suppose that SLE^ with ft < 4 has the restriction property. Then the distribution 
of all paths that are restricted not to hit D, and which are generated by SLE^ in H, is the same 
as the distribution of all paths generated by SLE„ in the domain H\D. SLE^ has the restriction 
property only for ft = 8/3 and no other values of ft [263]. We also expect such a property to 
hold for the continuum limit of SAWs, assuming it exists. 

The winding angle between the two endpoints of a finite 2D-SAW and indeed, a broader 
class of critical interfaces was first studied by Duplantier and Saleur |264| . Using conformal 
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invariance and nonrigorous Coulomb gas methods, they found that the distribution of winding 
angle approaches a Gaussian and they explicitly computed the winding variance as ~ ( 8 /gr) InL, 
where L is the distance between the end points of the walk, and g is a model dependent parameter 
which is 3/2 for SAW. The winding angle at a single endpoint relative to the global average 
direction of the curve is a Gaussian with variance ~ (4/g)lnL | 264| . It has been also found 
numerically | 265 | that the variance in the winding at typical random points along the curve was 
only 1/4 as large as the variance in the winding at the end points ~ (l/<?) luL. For SLE„ curves, 
the variance in the winding angle at the end point of the curve is shown to be kIuL |82| (the 
relation between the Coulomb gas parameter g and k can thus be given hy k = 4/g). 


4-2. Scaling limit and conformal invariance of percolation 

The existence of the conformally invariant scaling limit of percolation was first conjectured 
in based on experimental observations. This was then supported by some mathematical 
evidence provided for a different but related model, Voronoi percolation, which was proven |266| 
to be invariant with respect to a conformal change of metric. Using nonrigorous methods of 
CFT, Cardy could derive an exact limiting formula (21 1 for the crossing probability in a unit 
disc. Carleson has made an essential observation that this formula takes a particularly simple 
form when the domain is an equilateral triangle. In particular, for a percolation defined in an 
equilateral triangle A of side length 1 and vertices zi, Z 2 and Z 3 , and if 24 is on (23,21) at 
distance x € ( 0 , 1) from 23, then the crossing probability is simply Ps ((21,22), (23,24); A) = x. 
However, for years mathematicians were unable to rigorously justify the Cardy’s formula. 

In 2001, Smirnov [SUI proved that for site percolation on the triangular lattice, the limiting 
crossing probability exists which is conformally invariant and satisfies Cardy’s formula. Confor¬ 
mal invariance of the limit means that if the domain H, on which the percolation is dehned, to 
be conformally mapped onto any other domain H', such that 21 is mapped to 2j, 22 to 22, 23 to 
23 and 24 to Z 4 , then Ps ((21,22), (23, 24); H) = Ps((2(, 22), (23, 24); H'). The proof is based on 
the discovery of discrete harmonic functions which encode the crossing probability and converge 
to conformal invariants of the domain, in the scaling limit. 

Although we have been able to define the notion of a limiting crossing probability (though 
it may not exist at all), it was not a priori clear until about 2000 |49l 126711268| how to dehne 
a limiting percolation configuration i.e., a construction which does not involve limits of discrete 
systems. In fact, such a construction has been proposed by SLE„ curves as universal candidates 
for the scaling limits of macroscopic interfaces in 2D critical models. In order to show the 
existence and conformal invariance of the scaling limit for percolation, one can employ either its 
locality property or Cardy’s formula for crossing probabilities to show that k = 6 . Based on this 
observation Schramm conjectured [SS] that if percolation interface has a conformally invariant 
scaling limit, it must converge to SLEe—see also Fig. |12[ Smirnov has then outlined a proof 
for the conformal invariance of the full percolation conhguration. For more technical details of 
proof we refer to the original paper by Smirnov |80| or fFM\Tm\Tr 2 \. 

As an application of this convergence result, one can prove that the critical exponents for 
2D percolation exist, and their exact values can be computed, except for a and r, which are 
nevertheless listed here for completeness: 
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4-3. Percolation and magnetic models 

The analogy between bond percolation and conventional critical behavior in spin systems 
such as g—state Potts model, can be developed through a well-known mapping first discovered 
by Fortuin and Kasteleyn (FK) |27l 1281110311104| . The q—state Potts model is a generalization 
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Figure 12: Site percolation model on a square lattice in M. The fixed boundary condition is enforced 
at the lower boundary such that on the negative real half-line all the sites are unoccupied, while on the 
other half-line the sites are occupied. This imposes an interface (the perimeter in red) at the boundary 
of the spanning cluster (yellow colored) starting from the origin and ending at the upper boundary (the 
interface is defined uniquely by using the turn-right tie-breaking rule [269]). The external perimeter is 
shown by the dark solid line. In the scaling limit, as the lattice constant goes to 0, the perimeter and 
external perimeter then converge in distribution to SLEg and SLEg/ 3 , respectively. 


of the Ising model in which the spins at each site of a lattice can assume q possible spin values. 
It is defined by the lattice Hamiltonian H = ^^i(r), 3 (r'), where the sum is over nearest 

neighbors and Sij is the Kronecker delta. In the ferromagnetic case J > 0 (we set J = 1), 
the state in which all sites have the same spins minimize the energy and the system exhibits 
spontaneous magnetization at sufficiently low temperatures. There exists a critical temperature 
Tc at which the system undergoes a phase transition to the disordered phase. In two dimensions 
for q < 4 the phase transition at Tc is continuous. The partition function as a function of 
the inverse temperature /? is Z{j3) = Trexp 53(r),s(r'))! which apart from an overall 

unimportant constant may be rewritten as .E = Tr ((1 —p) +pSs(r),s(r')) with p = 1 — e~^ . 
Now every term in this expression is associated with a bond configuration in which there exists 
a bond for each term oc p and there is no bond for each term oc (1 — p). Sites connected by 
occupied bonds form clusters, and the Kronecker deltas force all the spins in each cluster to be 
in the same state. When the trace is taken over the spins, each cluster will have only one free 
spin and thus will give a factor q. Therefore one can write the partition function as a sum over 
conhgurations C of occupied bonds [273] 

Z = - p)^~^'’q^‘^ = percolation, (26) 

C 

where N is the total number of bonds in the lattice, Nb is the number of occupied bonds and 
Nc is the number of distinct clusters in C. This is the random cluster (or FK) representation 
of the g—state Potts model. Clearly the limit g —>■ 1 reproduces percolation with nontrivial 
correlations. 

From such a correspondence ( |26| ), one can build a dictionary which relates the thermody¬ 
namic quantities to the geometrical properties. For example, for the Ising model with q = 2, 
the spin-spin correlator is equal to the pair connectedness function of FK clusters i.e., the prob¬ 
ability that two sites of a distance r belong to the same FK cluster. The linear dimension of 
FK clusters diverges as the Ising correlation length, and the average (FK) cluster size diverges 
as the Ising susceptibility. Moreover, the average number of occupied bonds is proportional to 
the internal energy, and its fluctuations diverge as the Ising specific heat |274| . The presence of 
a spontaneous magnetization at T < Tc reflects the appearance of an infinite cluster at p > pc- 
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Figure 13: Illustration of GS clusters (left) versus FK clusters (right) for an Ising spin configuration. 
Different distinct clusters are shown in different colors. A GS cluster in a spin configuration, is a 
set of nearest-neighbor sites of like states. An FK cluster can then be constructed from a GS cluster 
by randomly assigning a bond between each pair of spins with a temperature dependent probability 
p := 1 — . The critical point on the square lattice is at /3c = ^ ln(l -h \/2) and pc = V2/{l + V2). 


For an Ising model on a 2D lattice with spin variables ct; = ±1, an alternative description of 
the partition function can be given by 2 = Trexp(—/35) with the action 5 = (1 — 

which only receives contributions from links across which the neighboring spins are anti-aligned. 
This means that if we have two adjacent clusters of opposite spins, the contribution to the action 


is proportional to the length of the boundary of geometric spin (GS) clusters—see Fig. 13 Thus 


the action can be considered as a sum over the self-intersecting connected admissible curves { 7 } 
ESH on the dual lattice weighted by their lengths L['y\ |275 |. l.e., 


Z = exp{ ^(-1) 
{7} 


"•(7)g-2/3i[7] j 


(27) 


where n{'y) cou nts the number of intersections of the immersed curve 7 . The topological term 
( — !)”('’') in (271 is essential since it actually distinguishes between the different intrinsic topolo¬ 


gies corresponding to a given extrinsic geometry and introduces cancellations between them to 
avoid over-counting of configurations. From here it is not so difficult to see that the continuum 
limit of this theory is a theory of free Majorana fermions that at the critical temperature Tc, 
becomes a conformal field theory with central charge c = 1/2 More recently the existence 
of conformal invariance and scaling limit of the 2D Ising model was rigorously proved. It has 
been shown that the geometric spin (GS) cluster interfaces as well as FK cluster interfaces in 
a 2D Ising model strongly converge to the SLE 3 and SLE 16 / 3 , respectively, in the scaling limit 
[ 851 1 ^ 1 ^ 8 ] . 


4-4- Dimensional reduction in criticality of a 3D Ising model 

The 3D Ising model has, so far, resisted an exact solution and may not be even computation¬ 
ally tractable |108| . Nevertheless much is known about its critical behavior, both analytically 
and numerically. 

Inspired by the formulation of the 2D Ising model as a theory of immersed curves, one may 
think of a possible extension of the same idea to the 3D model. There has been a lot of effort 
in the past [l03 EOl HHl HU EZSl ElOl EHJ |282l |283l |28l EMI EM] to reformulate the 3D 
Ising model in order to recast it as a string theory, i.e., as a theory of fluctuating membranes 
immersed in three dimensions. These attempts have been however stymied due to the difficulty 
in taking the continuum limit of formal sums over lattice surfaces. Part of the difficulty is that 


®The central charge c is the CFT parameter which is related to the SLE parameter k through c(k) = 
—— 12761 . Central charge is invariant under the duality k, —> 16/k. 
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in three dimensions the topological term (—(Z(E) is the nnmber of links where the closed 
lattice surface E intersects itself) oscillates very rapidly on the length scale of the lattice spacing. 

In this subsection we present a rather phenomenological approach based on application of 
percolation theory to give some evidence that the critical properties of a 3D Ising model are 
encoded in certain observables in a 2D cross section of the model. This suggests that is possible 
to employ the well-developed theory of immersed curves to study the Ising model in three 
dimensions. As discussed in the last subsection, the FK clusters in 5 —state Potts model always 
percolate right at the critical temperature Tc. The GS clusters in turn, do percolate at the 
same temperature Tc only in two dimensions. For example, in a 3D Ising model the FK clusters 
percolate exactly at the Curie point T^, while the percolation transition of the GS clusters, in 
contrast, occurs at some temperature Tp well below Tc |287 | and thus the 3D GS clusters do not 
capture the critical properties of the model. The results of extensive Monte Carlo study of 3D 
Ising model on a cubic lattice shows m that if one looks at an arbitrary 2D cross-section of the 
model, the GS clusters exhibit a percolation transition at a threshold which coincides exactly 
with the Curie point. It is also found numerically that the perimeter and the external perimeter 
of a GS cluster in a 2D cross-section at the Curie point satisfy the duality relation and their 
fractal dimensions and winding angle statistics are compatible, in the scaling limit, with SLE„ 
with ft = 5 and 16/5 respectively. This latter is in the same universality class as interfaces in 
tricritical Ising model in two dimensions. This numerical evidence may however pave the way 
to build a theoretical framework to understand the critical properties of the 3D Ising model. 

This observation on dimensional reduction of the criticality becomes more interesting if it 
would be a general feature of magnetic models independent of the dimensionality and the type 
of microscopic interactions. In fact some recent primary results |288| show that even for an Ising 
model in four dimensions, the GS clusters in a 2D cross-section of the model percolate exactly 
at the critical point of the original 4D model. 


5 . Percolation description of landscapes 

Percolation theory has been extensively applied to describe the properties of both artificial 
and natural landscapes. Percolation properties of the correlated surfaces as a model of wide 
range of artificial landscapes have been discussed in subsection |3.3| Our main focus in this 
section is mostly devoted to the statistical properties of natural landscapes and their modeling. 

The power spectrum S of linear transects of Earth’s topographjQ have a remarkable char¬ 
acteristic scaling relation with the wave number k as S(k) ~ k~^'^ with the exponent /3c = 2, 
over a wide range of scales |289l 12901129111292| . Similar scaling relations have been identified 
in Earth’s bathymetry (i.e., the underwater equivalent to topography) |293j . the topography of 
natural rock surfaces |294| . and the topography of Venus |295| . Such a power-law spectrum in 
the topography is responsible for the appearance of various self-similar patterns on Earth, e.g., 
fractal coastlines [296], the radiation fields of volcanoes [29711298] . crustal density and gravity 
|299| . geomagnetism |30U| . and surface hydrology such as in the river basin geomorphology | 3U1 |. 
Although environmental parameters such as erosion seem to play an important role in shaping 
the coastlines, drainage basins and watersheds, the observation of scale-invariant topography on 
Venus, however, indicates that fractal topography can be formed without erosion. 

The exponent f3c is related to the Hurst exponent H in fractional Brownian motion (fBm) via 
Pc = 2JT -I- 1, thus suggesting H ~ 0.5 for Earth’s topography. However, further surveys based 


^The power spectrum S{k) is defined as the square of the coefficients in a Fourier series representation 
of the transect, which measures the average variation of the function at different wavelengths. For totally 
uncorrelated adjacent data points S{k) is a constant, while for strongly correlated ones relative to points 
far apart, it will be large at small k (long wavelengths) and small at large k (short wavelengths). 
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on the fBm model |29U| of topography or bathymetry revealed a more complex multifractal 
structure of Earth’s morphology giving rise to distinct scaling properties of oceans, continents, 
and continental margins describes by = 0.46,0.66 and 0.77, respectively | 302| . 

The other characteristic feature of Earth’s topography is its bimodal distribution |303| which 
reflects the topographic dichotomy of continents and ocean basins. It has a clear discrepancy 
with Gaussian models of topography, also consistent with the Mandelbrot’s observation | 304 | . 
The positive correlation between elevation and slope seen on Earth (i.e., the steepness increases 
with the height) is not also predicted by a model with a Gaussian distribution, implying that 
the global topography of Earth is not easily amenable to modeling. 

Nevertheless, percolation theory has been recently applied to describe the global topography 
of Earth [21] , in which the critical point indicates the present mean sea level. Moreover, different 
models based on percolation theory have been proposed to describe the statistical properties 
of regional features on Earth such as coastlines river basins and drainage networks 

[3061 13071 13081 13091 13101 1311| . and watersheds | 1151 11161 11291 I312| . Percolation theory has 
been also successful to help understand other phenomena on Earth. It has been demonstrated 
ESI that sea ice exhibits a percolation transition at a critical temperature above which brine 
carrying heat and nutrients can move through the ice, whereas for colder temperatures the ice 
is impermeable. Percolation also serves as an attractive mechanism to explain core formation 
in Earth |314ll3l5| . 

5.1. Fractal geometry of coastlines 

Coastlines were among the first natural systems that have been quantitatively characterized 
when Mandelbrot computationally analyzed their fractal geometry |296| . In fact the geometrical 
irregularity of the coastlines helps damping the sea waves and decreasing the average wave 
amplitude. Affected by the sea eroding power, an irregular morphology evolves at the rocky 
coast until a self-stabilization with the wave amplitude is established. A simple model of such 
stabilization has been studied |20| in which the fractal geometry of the coastline plays the 
role of a morphological attractor directly related to percolation geometry. Dynamics of the 
model spontaneously leads to a stationary fractal geometry with a dimension very close to 
4/3 independent of the initial morphology, in agreement with that is observed on real coasts 
[29611316] . This fractal dimension is also consistent with that of the external perimeter of the 
spanning cluster in a 2D critical percolation. Two general erosion mechanisms are considered in 
the 2D model, i.e., a rapid mechanical erosion and a slow chemical weakening. It has been shown 
that when the model involves both processes, a dynamic equilibrium is reached that changes 
the shape of the coast but preserves its fractal properties. 

The effect of spatial long-range correlations in the lithology of coastal landscapes on the 
fractal properties of the coastlines has then been addressed in [305] . In fact, due to the endogenic 
processes like volcanic activity, earthquakes, and tectonic processes originating within Earth that 
are mainly responsible for the very long-wavelength topography of Earth’s surface, one naturally 
expects that lithological properties of coastal landscapes would be in general heterogeneous as 
well as long-range correlated in space. Moreover, a multitude of fractal dimensions has been 
measured for real coastlines of different landscapes m\- Thus self-similar geometry of coastlines 
should emerge from an intricate interplay between these landscape properties and the sea force. 
The results of a simple invasion model [305] indicates that a critical sea force fc exists at 
which the coastline exhibits self-similarity with fractal dimension depending on Hurst exponent. 
The dominant 4/3 fractal dimension was obtained for uncorrelated landscapes. For f < fc the 
coastline is rough but not fractal and the eroding process stops after some time, while for f > fc, 
erosion is perpetual leading to a self-affine coastline which belongs to the Kardar-Parisi-Zhang 
(KPZ) [318] universality class. 

As discussed in Sec. [^ the external perimeter of critical percolation clusters with fractal 
dimension 4/3 are proven to have a conformally invariant limit described by SLEg/a. Some 
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numerical evidence of such a strong symmetry has been also reported for rocky coastlines with 
fractal dimension 4/3 | 1U2 |. These coastlines are therefore shown to be statistically equivalent 
to the external perimeter of percolation clusters or that of planar random walk. The conformal 
invariance can then be used to predict the statistics of the flux of pollutants diffusing over 
shorelines. This flux has been characterized by a strongly intermittent spatial distribution 
which can vary dramatically between locations just a few hundred meters apart. 

Strong evidence of conformal invariance property has been also presented for the iso-height 
lines (like coastlines) of both artificial landscapes in the KPZ universality class |93l 1951 I97| 
and experimentally grown surfaces |92]. In the KPZ universality class, the iso-height lines are 
characterized by a fractal dimension of 4/3 with the same conformal invariant properties as the 
external perimeter of critical percolation clusters, compatible with SLEg /3 curves in the scaling 
limit. Such an analogy may lead to an alternative description of the coastlines. 

5.2. Statistical properties of watersheds 

The watershed is defined as the line which separates adjacent drainage basins—see Fig. 
Based on observation of natural watersheds, it has been claimed that they should have a fractal 
structure The self-similarity of watersheds has then been justified numerically for both 

natural and artificial landscapes [TT51 [Uni 1320]. Watersheds have been shown to be related 
to a family of curves appearing in different contexts, e.g., bridge percolation | 192| . polymers 
in strongly disordered media ism, optimal path cracks |322| . and fracturing process |323| . To 
determine the watershed lines on real or artificial landscapes which are usually in the form of 
Digital Elevation Maps (DEM), consisting of discretized elevation fields, one can use an iterative 
application of an invasion percolation procedure 

The fractal dimensions of watershed lines in 2D and 3D were estimated ESI to be d/ = 
1.2168 ± 0.0005 and 2.487 ± 0.003, respectively, for uncorrelated artificial landscapes. In two 
dimensions, however, the measured fractal dimensions for natural landscapes obtained from data 
provided by satellite imagery |324 |. fall into the range 1.10 < df < 1.15. This may imply the 
necessity of considering spatial correlations in computations. When the long-range correlations 
characterized by the Hurst exponent H were introduced |320| . a monotonic decrease of d/ with H 
has been observed, and the agreement with the observation achieved for 0.3 < H < 0.5 (although 
this range of H seems to be out of that is observed for continents and continental margins 
[302]). Moreover, it has been shown that small and localized perturbations like landslides 
or tectonic activities, can have a large and non-local impact on the shape of watersheds. It 
is also discussed in | 312| that the fractal dimension obtained in 2D for uncorrelated artificial 
landscapes is intriguingly close to the fractal dimension of the largest cluster boundary in two 
models of explosive percolation on a lattice, i.e., the largest cluster and Gaussian models HU. 

Watersheds are shown |99| to be among the rare examples of physical systems described by 
SLEk curves with k < 2. It has been numerically shown that, in the scaling limit, the watershed 
line exhibits conformally invariant properties compatible with SLE„ with k = 1.734 ± 0.005. 

5.3. The present mean-sea level on Earth 

The ubiquitous scale invariant features on Earth have endowed theoretical interest on the 
assumption that they may reveal prevalence of some underlying feature [3251132611327| . This is 
still an open question if there exists a clear relationship between the quantitative properties of 
landscapes and the dominant geomorphologic processes that originate them. Although such a 
relationship is established for some of regional features, the global topography in comparison, 
has received less attention. 

The appearance of scale and conformal invariance property in statistical models like perco¬ 
lation is a specific feature of criticality. This can be regarded as a motivating issue to search for 
an underlying mechanism that possibly explains the emergence of fractal geometries on various 

36 




landscapes. For instance, it has been shown for an ensemble of experimentally grown surfaces 
|34 | that there exists a critical level height at which a percolation transition occurs. This may 
elucidate the earlier observation [92] of conformal invariant iso-height lines on these samples. 



Figure 14: Schematic illustration of the continental aggregation by decreasing the sea level from top 
{h = -1-100 m) to bottom {h = —80 m). This shows a remarkable percolation transition at the present 
mean sea level around which the major parts of the landmass join together. 

A percolation description of the global topography of Earth is recently presented m in 
which a dynamic geoid-like level is defined as an equipotential spherical surface as a counterpart 
of the percolation parameter. When the hypothetical water level is decreased from the highest 
to lowest available heights on Earth, there occurs a geometrical phase transition at a certain 
critical level he around which the most parts of landmass join together—see Fig. |14| The most 
remarkable observation is that the critical level he coincides with the present mean sea level 
h — 0 on Earth. The criticality of the current sea level justifies the appearance of the scale (and 
conformal) invariant features on Earth. This may also uncover the important role that is played 
by water on Earth and shed new light on the tectonic plate motion. 



Figure 15: Map of tectonic plates compared with different disjoint islands for the sea level at h = -1-100 m. 
Every disjoint landmass is approximately surrounded by a major plate boundary. 

According to the plate tectonic theory, the outer portion of Earth is made up of a number 
of distinct plates (Fig. |15[) which move relative to each other. This motion is responsible for the 
major topographical features such as creation of oceans and pushing up mountain ranges. The 
open question motivated by this work is whether such an observed criticality plays the role of a 
geometric attractor for tectonic motions through geological time. 
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6 . Conclusions 


Through the current review we have outlined some basic properties and recent advances 
of percolation theory as well as some of its recent applications which can be summed up as 
the following statements. Percolation theory and its applications span a wide area of science 
ranging from social and network sciences to string theory and particle physics as well as various 
branches of probability theory in mathematics. The percolation models are mostly governed by 
very simple rules yet with a fascinating mathematical structure and fundamental features. It has 
an uncorrelated microscopic structure, but a long-ranged correlated geometry can emerge which 
governs the system. All its properties can be described by a number of geometric observables 
with some characterizing features distinguishing between subcritical, critical and supercritical 
phases. The criticality has a rich fractal structure and remarkable underlying scaling laws 
defined by some universal critical exponents. Percolation theory simultaneously benefits from 
exact conjectures raised by physical insights on the one hand, and rigorous mathematical proofs 
on the other hand. It has a very robust nature against small perturbations but ready to play a 
role in a completely different scene under sufficiently large modihcations. With no interactions, 
it may be discovered in the heart of strongly interacting systems While it is applied to 
formulate the critical behavior of a system in terms of some appropriate observables, it can be 
viewed at the same time as another theory in a lower dimension describing the same system 
with appropriately well defined observables |31| . Many exact results have been obtained but 
there still exist many open challenges in the field |40| . 
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